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Abstract 

This  thesis  presents  an  airborne  radar  model  for  non-uniformly  spaced  antenna 
arrays.  An  airborne  radar  model  for  non-uniformly  spaced  arrays  is  important  because 
the  performance  of  arbitrary  array  configurations  can  be  analyzed.  Non-uniformly 
spaced  arrays  could  potentially  possess  advantages  which  are  unobtainable  with  uni¬ 
formly  spaced  arrays.  First,  an  element  position  matrix  is  defined  to  accommodate 
arbitrary  element  locations.  The  new  element  position  matrix  produces  changes  to 
the  time/phase  delay  to  each  element,  the  spatial  frequencies,  the  steering  vectors, 
the  space-time  snapshots,  and  the  covariance  matrices.  The  joint  domain  localized 
(JDL)  and  factored  time-space  (FTS)  STAP  algorithms  are  also  updated  so  they 
are  compatible  with  the  model  for  non-uniformly  spaced  arrays.  The  non-uniformly 
spaced  array  radar  model  presented  is  verified  with  previous  models  presented  by 
Jaffer  [6],  Ward  [13],  and  Hale  [5].  Past  research  involving  non-uniform  arrays  mostly 
addressed  circular  arrays,  particularly  the  ultra  high  frequency  (UHF)  Electronically 
Scanned  Array  developed  by  the  Office  of  Naval  Research.  Therefore,  the  model  for 
non-uniform  arrays  is  used  to  evaluate  performance  of  two  different  circular  arrays, 
one  containing  24  elements  and  one  containing  15  elements.  The  circular  array  per¬ 
formance  is  compared  to  that  of  a  6  x  6  uniformly  spaced  planar  array  having  the 
same  dimensions.  Signal-to-interference  plus  noise  ratio  (SINR)  loss  and  detection 
probability  metrics  are  used  to  evaluate  array  performance.  Array  failure  is  also 
examined  using  output  SINR  plots  for  each  array.  Simulation  results  indicate  that 
non-uniformly  spaced  arrays  such  as  the  circular  array  can  potentially  be  used  in  air¬ 
borne  radar  systems  in  place  of  uniformly  spaced  arrays.  As  a  result,  array  costs  may 
be  reduced  by  using  fewer  elements  without  sacrificing  detection  capabilities. 
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An  Airborne  Radar  Model  For 
Non-Uniformly  Spaced  Antenna  Arrays 

I.  Introduction 

Non-uniformly  spaced  antenna  arrays  are  studied  because  they  may  be  used  with 
space-time  adaptive  processing  (STAP)  techniques  in  airborne  radar  systems  in  the 
near  future.  Non-uniformly  spaced  arrays  are  a  topic  of  interest  since  they  may  have 
certain  advantageous  characteristics  not  possessed  by  uniform  arrays.  For  example, 
circular  arrays  can  be  electronically  scanned  360°  in  azimuth,  whereas  uniform  arrays 
must  be  mechanically  steered.  Circular  arrays  are  becoming  more  popular  since  they 
can  be  easily  implemented  in  the  nose  of  an  aircraft  [7] .  Some  non-uniform  array  con¬ 
figurations  may  also  provide  lower  sidelobes  than  uniformly  spaced  arrays.  Another 
possible  advantage  of  non-uniformly  spaced  arrays  is  the  use  of  fewer  elements.  Array 
costs  can  be  reduced  if  fewer  elements  are  arranged  non-uniformly  while  maintaining 
the  same  detection  performance  as  uniformly  spaced  arrays.  In  some  cases,  uniform 
spacing  of  elements  may  not  be  feasible.  There  might  be  limited  space  on  a  platform, 
so  an  antenna  array  may  have  to  conform  around  certain  components  or  other  obsta¬ 
cles  on  the  platform.  Non-uniformly  spaced  arrays  are  also  investigated  to  see  if  they 
are  more  or  less  susceptible  to  element  failure  than  the  uniformly  spaced  array. 

1.1  Objectives 

The  primary  objective  of  this  thesis  is  to  develop  a  radar  model  for  arrays  with 
any  element  spacing,  uniform  or  non-uniform.  The  radar  data  model  for  linear  arrays 
developed  by  Ward  [13]  and  Jaffer  [6]  and  extended  to  planar  arrays  by  Hale  [5]  is  mod¬ 
ified  to  include  non-uniform  element  spacing  for  linear,  planar,  and  three-dimensional 
arrays.  However,  the  model  derived  in  this  thesis  does  not  account  for  shadowing 
effects  that  may  occur  in  the  three-dimensional  array.  Shadowing  occurs  when  some 
elements  in  the  array  do  not  receive  the  incoming  radar  wave  due  to  their  location. 
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The  modifications  to  the  radar  model  begin  by  redefining  the  element  position  vector 
as  an  element  position  matrix.  The  new  element  position  matrix  trickles  through  the 
radar  model,  bringing  about  many  changes.  The  new  element  position  matrix  affects 
the  time/phase  delay  to  each  element,  the  spatial  frequencies,  the  steering  vectors, 
the  space-time  snapshots,  the  array  pattern,  and  the  covariance  matrices.  The  joint 
domain  localized  (JDL)  and  factored  time-space  (FTS)  STAP  algorithms  are  also 
modified  to  reflect  the  changes  due  to  non-uniform  element  spacing.  These  changes 
affect  output  signal-to-interference  plus  noise  ratio  (SINR),  SINR  loss,  and  detec¬ 
tion  probability.  Grating  lobes  may  also  result  from  non-uniform  element  spacing, 
depending  on  the  array  configuration. 

Once  a  model  for  non-uniform  element  spacing  is  developed,  the  next  objective 
is  to  ensure  it  is  correct.  To  verify  it  is  correct,  uniform  element  spacing  is  inserted 
into  the  model  and  compared  to  the  radar  model  of  [5,6,13].  The  two  models  are 
mathematically  equivalent  when  the  uniformly  spaced  array  is  linear.  The  two  mod¬ 
els  are  structured  differently  for  the  uniformly  spaced  planar  array  case.  Therefore, 
simulations  are  performed  to  verify  the  same  results  are  produced  from  each  model. 

1.2  Performance  Evaluation 

This  thesis  also  intends  to  determine  the  impact  of  non-uniform  element  spacing 
on  adaptive  interference  suppression  techniques.  The  performance  of  non-uniformly 
spaced  arrays  and  uniformly  spaced  arrays  is  compared  to  determine  if  certain  non- 
uniform  array  configurations  can  be  used  in  place  of  uniform  arrays.  Two  different 
non-uniform  arrays  are  used  for  experimentation  to  characterize  their  performance. 
The  non-uniform  arrays  used  include  a  15  element  circular  array  and  a  24  element 
circular  array.  The  STAP  techniques  used  in  the  simulations  are  the  JDL  and  FTS 
methods.  JDL  and  FTS  are  used  because  they  are  practical  STAP  methods  that 
can  be  implemented  in  real-world  radar  systems  due  to  their  reduced  sample  support. 
The  matched  filter  is  also  used  in  simulations  as  the  ideal  performance  bound  for  each 


2 


STAP  method.  The  performance  of  each  array  is  compared  to  evaluate  the  impact  of 
different  element  configurations. 

Array  performance  is  evaluated  based  on  a  variety  of  simulations.  The  antenna 
pattern  is  plotted  first  to  examine  the  energy  radiated  from  the  array.  The  antenna 
pattern  shows  features  such  as  sidelobes,  backlobes,  and  nulls.  Next,  the  clutter  co- 
variance  matrix  eigenvalues  are  plotted  to  examine  the  element  configuration  impact 
on  clutter  rank.  Covariance  matrices  describe  the  expected  signal  power  returns  from 
sources  such  as  clutter,  noise,  jammers,  or  targets.  A  matrix  with  high  rank  means 
there  is  weak  correlation  between  samples.  For  example,  the  white  noise  covariance 
matrix  is  full  rank  because  white  noise  is  completely  uncorrelated  from  sample  to 
sample.  The  jammer  and  clutter  matrices  have  low  ranks  due  to  correlation  between 
samples.  The  physical  meaning  of  low  clutter  matrix  rank  is  new  incoming  data 
samples  do  not  offer  much  more  information  about  the  clutter  environment  than  pre¬ 
vious  samples.  The  power  spectra  of  each  array  is  plotted  using  the  signal  match  and 
minimum  variance  estimator  (MVE)  algorithms  to  examine  the  interference  environ¬ 
ment.  The  signal  match  method  shows  the  interference  as  illuminated  by  the  array. 
Therefore,  interference  due  to  sidelobes  and  grating  lobes  is  contained  in  the  signal 
match  power  spectra.  The  MVE  power  spectra  simply  shows  where  the  interference 
is  physically  located  relative  to  the  radar.  The  performance  of  each  array  using  STAP 
techniques  is  finally  evaluated  with  SINR  loss  and  detection  probability  plots.  SINR 
loss  is  the  actual  output  SINR  referenced  to  the  maximum  output  signal-to-noise  ra¬ 
tio  (SNR).  SINR  loss  is  plotted  to  compare  array  performance  since  each  array  has 
a  different  number  of  elements,  resulting  in  a  different  output  SINR.  SINR  loss  is 
plotted  versus  normalized  Doppler  frequency  for  the  matched  filter,  FTS,  and  JDL 
methods.  Detection  probability  is  plotted  versus  input  SINR  per  element  per  pulse 
for  the  adaptive  matched  filter  (AMF),  JDL,  and  FTS  methods.  The  SINR  loss  and 
detection  probability  plots  characterize  the  performance  of  each  array  for  each  STAP 
method.  Finally,  array  failure  is  briefly  examined  using  output  SINR  plots  for  each 
array  with  0%,  10%,  and  20%  element  failure  rates. 
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1 . 3  Organization 

This  thesis  is  broken  down  into  six  chapters.  Chapter  I  introduces  the  problem 
and  research  goals.  A  brief  overview  of  how  the  problem  is  to  be  solved  is  also 
discussed.  Chapter  If  discusses  the  radar  data  models  presented  in  [5,6,13].  These 
models  are  for  uniformly  spaced  arrays.  A  few  different  STAP  algorithms  discussed 
in  [5,  12,  13]  are  also  mentioned  in  Chapter  II  since  they  are  used  in  simulations 
in  the  following  chapters.  Chapter  III  shows  the  development  incorporating  non- 
uniform  element  spacing  into  the  radar  model  of  [5].  Spatial  frequencies,  steering 
vectors,  space-time  snapshots,  covariance  matrices,  and  two  STAP  algorithms  are 
updated  to  reflect  the  necessary  changes.  The  model  presented  in  Chapter  III  is 
verified  in  Chapter  IV  to  ensure  it  is  correct.  Verification  is  done  by  inserting  uniform 
element  spacing  into  the  model  presented  in  Chapter  III  and  comparing  it  to  the  model 
in  [5,6,13].  First,  the  uniform  linear  array  is  verified  mathematically.  Then,  the  planar 
array  is  verified  via  simulation  since  the  variables  are  structured  differently  between 
each  model.  Simulation  plots  and  results  are  presented  in  Chapter  V  to  evaluate  the 
performance  of  each  array.  SINR  loss,  output  SINR,  detection  probability,  clutter 
rank,  and  power  spectra  plots  are  the  simulations  used  for  evaluation.  Chapter  VI 
draws  conclusions  based  on  the  simulation  results.  Future  research  is  also  discussed 
in  this  chapter. 

1.4  Notation 

All  notation  in  this  thesis  follows  the  standard  technical  writing  convention. 
Scalar  values  are  symbols  in  regular  font.  For  example,  6  is  the  elevation  angle  and 
cj)  is  the  azimuth  angle.  Column  vectors  are  represented  by  lower  case  bold  symbols. 
The  temporal  steering  vector  b  and  the  space-time  steering  vector  v  are  examples  of 
column  vectors.  Matrices  are  represented  by  upper  case  bold  symbols.  For  example, 
R  is  the  covariance  matrix.  Identity  matrices  always  have  subscripts  denoting  their 
dimensions.  An  identity  matrix  of  size  NM  x  NM  is  written  Inm-  The  Hermitian 
is  the  complex  conjugate  transpose  and  is  denoted  by  a  superscript  H  as  in  HH .  A 
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transformed  vector  or  matrix  has  a  tilde  over  it.  The  vector  x  is  X  projected  into 
a  different  space.  Estimated  quantities  and  unit  vectors  have  hats  over  them.  For 
example,  R  is  an  estimate  of  the  covariance  matrix  R  while  a:  is  a  unit  vector  in 
the  x  direction.  Careful  attention  must  be  paid  to  the  context  to  distinguish  which 
operation  the  hat  is  being  used  for.  The  expected  value  operator  is  written  as  £{•}. 
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II.  Past  Research  In  STAP  For  Airborne  Radar 


This  chapter  examines  past  research  efforts  involving  the  model  for  an  airborne 
radar  and  STAP  algorithms  for  that  model.  A  model  for  radar  system  analysis 
is  presented  in  [6, 13]  for  a  sidelooking  uniform  linear  array.  That  model  was  later 
extended  in  [5]  for  a  uniform  planar  array  containing  P  rows  of  uniformly  spaced  N 
element  linear  arrays.  That  radar  model  is  discussed  here  because  it  is  modified  in 
the  following  chapter  to  allow  non-uniform  element  spacing. 

2.1  Airborne  Radar  Model 

2.1.1  Array  Geometry.  The  radar  model  in  [5]  begins  with  the  physical 
geometry  of  the  antenna  shown  in  Fig.  2.1.  The  elements  are  divided  into  P  rows 
along  the  2-axis  and  N  columns  along  the  x-axis.  The  space  between  each  element  is 
denoted  as  dx  along  the  x-axis  and  dz  along  the  2-axis.  The  location  of  element  np  is 
defined  as 

d np  =  —ndxk  -  pdzz,  (2.1) 

where  x  and  z  are  unit  vectors  in  the  x  and  2  directions,  respectively.  The  array 
geometry  affects  every  aspect  of  the  radar  model. 


Figure  2.1:  Sidelooking  uniformly  spaced  planar  array  lying 

in  the  x2-plane.  The  reference  element  is  located  at  the  origin. 
The  aircraft  velocity  vector  is  along  the  positive  x-axis. 
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2.1.2  Received  Signal.  The  received  signal  time/phase  delay  to  each  array 
element  changes  based  on  the  array  geometry.  The  time  delay  to  the  npth  element  is 


'  np 


k(6>,  0)  •  dnp 
c 

—ndx  cos  6  sin  0  —  pdz  sin  9 
c 


(2.2) 


The  resulting  spatial  frequency  d  for  a  target  located  at  azimuth  angle  0  and  elevation 
angle  6  also  changes  based  on  the  array  geometry.  The  spatial  frequency  is  divided 
into  an  ^-component 


'0 


X 


dx  cos  6  sin  0 


A 


O 


(2.3) 


and  ^-component 


dz  = 


dz  sin  9 

A0 


(2.4) 


when  the  array  is  uniformly  spaced.  A  uniformly  spaced  planar  array  can  have  two 
spatial  steering  vectors,  one  for  azimuth  and  one  for  elevation,  to  locate  a  target  in 
space.  The  azimuth  steering  vector  a  and  elevation  steering  vector  e  are  defined  as 


a(tf*)  =  [1  e'2’""  ■  ■  ■  e-’2’('v~1),’-]T 


(2.5) 


and 

e(dz)  =  [1  ej2^z  ■  ■  ■  eV<P-RV*Y- 
The  temporal  steering  vector 


b(w)  =  [1  e22”2’  •  ■  ■  eiMM~  1)"]T 


(2.6) 


(2.7) 


is  used  with  the  spatial  steering  vectors  to  construct  the  3D  steering  vector  to  a 
specific  azimuth  angle,  elevation  angle,  and  normalized  Doppler  frequency  do  as 


v(dz,iv,  '9x)  =  e(dz)  <g)  b(<n)  (8)  a (dx). 


(2.8) 
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The  Kronecker  Product  ®  is  described  in  Appendix  A  and  M  is  the  number  of  pulses. 
The  space-time  snapshot  Xt  representing  the  target  return  is  simply  the  3D  steering 
vector  scaled  by  a  complex  value 

Xt  =  <xtv(‘&z,  u^x)-  (2.9) 

2.1.3  Array  Factor.  The  array  factor  is  necessary  to  construct  the  correct 
transmit  pattern  from  the  array.  The  array  factor  is  required  for  clutter  and  jammer 
modelling  since  it  dictates  the  energy  amplitude  and  location  in  azimuth  and  elevation. 
The  array  factor  is  the  sum  of  all  weights  over  the  array  and  is  defined  as 

W(e,(f>)  =  [e(tf2)  <g>  a(^)]Tl  (2.10) 

in  [5]  where  1  is  a  vector  of  ones.  Equation  (2.10)  is  the  array  factor  for  a  mainbeam 
steered  to  boresight.  The  mainbeam  can  be  steered  to  any  azimuth  and  elevation  by 

W(e,(j))  =  [e(d2)  ®  a(^)]H[e«)  a(^)],  (2.11) 

where  rt'x  and  '0'z  contain  the  mainbeam  azimuth  and  elevation  location.  See  [1,11]  for 
more  information  on  array  beamforming. 

2.1.4  Jammer  Model.  The  jammer  modelled  in  [5,6,13]  is  a  barrage  noise 
jammer.  This  barrage  noise  jammer  model  is  briefly  mentioned  here  and  is  discussed 
in  greater  detail  in  Chapter  III.  The  jammer  space-time  snapshot  is 

Xj  =  e(d2)  (8)  otj  <g>  a(dx)  (2.12) 


where 


otj  —  [ckq  ot  1  ■  •  •  olm— 1  T 
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(2.13) 


is  a  vector  representing  the  complex  random  jammer  amplitudes  received  from  each 
pulse.  A  barrage  noise  jammer  has  no  temporal  correlation  since  its  amplitude  varies 
from  pulse  to  pulse.  The  jammer  covariance  matrix  is  constructed  from  its  space-time 
snapshot  as 

R,  =  4XjX?  }■  (2.14) 

2.1.5  Clutter  Model.  The  clutter  model  in  [5,6,13]  is  briefly  mentioned  here 
since  it  is  discussed  in  greater  detail  in  the  following  chapter.  A  constant  gamma 
model  for  ground  clutter  is  the  clutter  model  basis.  The  clutter  space-time  snapshot 

is 

Nr- 1  Nc- 1 

Xc  =  51  a^e(^z)  ®  b(cAfc)  ®  a (rfx),  (2.15) 

i=0  k= 0 

where  cq*,  is  the  clutter  patch  amplitude  and  dz  and  dx  contain  the  elevation  angle 
9-i  and  azimuth  angle  dk  to  the  clutter  patch.  The  clutter  patch  normalized  Doppler 
frequency  Co  is 

Oik  =  /c(g;0fc)  =  UOiMTr 
Jr 

=  =  (2.16) 

dx 

where  fc  is  the  Doppler  frequency,  fr  is  the  pulse  repetition  frequency  (PRF),  Tr  is  the 
pulse  repetition  interval  (PRI),  and  va  is  the  aircraft  speed.  Equation  (2.16)  shows  the 
normalized  Doppler  frequency  dependence  on  azimuth  angle  with  a  slope  /3  describing 
the  number  of  half-interelement  spacings  dx  crossed  during  one  PRI.  The  value  of  [3 
determines  how  difficult  it  is  to  suppress  clutter.  If  f3  is  greater  than  one,  clutter 
spatial  frequencies  are  spread  over  a  frequency  band  greater  than  the  PRF,  resulting 
in  Doppler  ambiguous  clutter  returns.  In  other  words,  different  clutter  patches  can 
have  the  same  Doppler  frequency  as  each  other.  The  final  portion  of  the  clutter  model 
is  the  covariance  matrix 

=  UXoX?}.  (2.17) 
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The  clutter  model  is  important  because  it  describes  the  interference  from  the  envi¬ 
ronment  around  the  radar  platform. 


2.2  Performance  Metrics 

A  few  different  performance  metrics  are  used  to  characterize  the  performance  of 
an  array.  These  metrics  include  output  SINR,  SINR  loss,  and  detection  probability. 
The  signal  match  and  MVE  methods  are  two  different  spectral  estimation  techniques 
used  to  examine  the  interference  environment. 


2.2.1  Output  SINR.  Output  SINR  is  described  in  [13]  as  a  common  metric 
used  to  characterize  processor  performance.  It  is  a  ratio  of  output  target  power  to 
interference  plus  noise  power.  SINR  is  given  as 


SINR 


cr%.\wHv\2 

wffRw 


(2.18) 


in  [5],  where  is  the  target  SNR  per  element  per  pulse,  a2  is  the  noise  power  per 
element  per  pulse,  and  w  is  the  weight  vector.  SINR  is  plotted  as  a  function  of 
normalized  target  Doppler  frequency  since  target  velocity  is  not  known.  Output  SINR 
changes  based  on  the  number  of  elements  in  the  array.  Arrays  with  more  elements 
have  a  higher  output  SINR  than  arrays  with  fewer  elements. 


2.2.2  SINR  Loss.  SINR  loss  is  another  performance  metric  described  in 
[13].  It  is  valuable  because  it  shows  the  losses  due  to  interference  suppression,  filter 
straddling  losses,  and  taper  losses.  Therefore,  interference  suppression  performance 
of  arrays  with  differing  numbers  of  elements  can  be  compared.  SINR  loss  is  defined 
as 

SINR 
is,NR  “  SNR? 

where 

SNR0  =  £tMNP  (2.20) 
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is  the  maximum  output  SNR  limited  only  by  noise,  the  number  of  pulses,  and  the 
number  of  array  elements.  SNR0  is  the  optimum  SNR  achieved  without  interference. 


2.2.3  Signal  Match  Power  Spectra.  The  signal  match  power  estimator  is 
discussed  in  [7].  It  shows  the  interference  due  to  the  surrounding  environment  and  the 
antenna  pattern  without  adaptivity.  Interference  illuminated  by  the  antenna  pattern 
sidclobes  shows  up  in  the  signal  match  power  spectra.  The  normalized  output  power 
for  the  signal  match  method  is 


P  = 


vH  Rv 


(2.21) 


The  signal  match  method  is  a  Fourier  transform  when  element  spacing  is  uniform. 


2.2.4  Minimum  Variance  Estimator  Power  Spectra.  The  other  spectral 
estimation  technique  is  a  high  resolution  technique  called  the  MVE  [7].  The  MVE 
algorithm  decomposes  the  signal  into  single  peaks  and  shows  where  the  interference 
is  physically  located  relative  to  the  radar.  Amplitude  information  under  the  MVE  is 
meaningless  because  its  high  resolution  capabilities  introduce  instabilities,  causing  a 
loss  in  amplitude  estimation  ability.  Therefore,  MVE  simply  tells  you  if  a  signal  is 
present  or  not.  The  output  power  using  the  MVE  method  is 


P  =  (vhR-\)-\ 


(2.22) 


2.3  STAP  Methods 

A  few  different  STAP  methods  are  examined  for  interference  suppression  capa¬ 
bility.  The  matched  filter  and  AMF  attain  the  best  possible  performance,  meaning 
they  do  the  best  job  at  suppressing  interference.  However,  the  matched  filter  and 
AMF  are  not  practical  because  they  require  too  much  sample  support.  Reed  [9] 
determined  that  twice  the  degrees  of  freedom  (DOF)  sample  support  vectors  are  re¬ 
quired  to  obtain  performance  within  3  dB  of  optimal.  Each  practical  STAP  method 
uses  reduced  DOF  to  achieve  performance  comparable  to  the  matched  filter  and  AMF. 
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FTS  and  JDL  are  two  practical  STAP  methods  examined  here  since  they  are  used 
for  simulations  in  Chapters  IV  and  V.  The  JDL  and  FTS  algorithms  are  presented 
in  [5]  for  uniformly  spaced  planar  arrays  and  are  updated  in  Chapter  III  to  reflect  the 
changes  caused  by  non-uniform  element  spacing. 

2.3.1  Matched  Filter.  The  matched  filter,  also  called  the  optimum  filter, 
can  be  found  using  the  likelihood  ratio,  the  maximum  SINR,  or  the  Wiener  filter 
method  [2].  The  matched  filter  weights  are  found  by  projecting  the  space-time  steering 
vector  onto  the  covariance  matrix  inverse, 

w  =  R  V  (2.23) 

The  matched  Liter  is  considered  non-adaptive  since  the  covariance  matrix  R  of  incom¬ 
ing  space-time  snapshots  is  known  prior  to  receiving  the  incoming  data.  Detection 
probability  is  maximized  by  maximizing  SINR  using  the  matched  filter. 

2.3.2  Adaptive  Matched  Filter.  The  AMF  is  similar  to  the  matched  filter, 

except  the  covariance  matrix  of  incoming  space-time  snapshots  is  not  known  prior  to 

receiving  the  data.  Therefore,  the  covariance  matrix  is  estimated  from  the  data  al¬ 
ready  received  using  the  Maximum  Likelihood  Estimate  (MLE)  [9]  and  is  represented 
as  R.  The  AMF  weights  are 

w  =  R  v.  (2.24) 

Incoming  data  must  be  estimated  since  it  is  constantly  changing  due  to  aircraft  mo¬ 
tion,  heterogeneous  clutter,  jamming,  etc.  The  AMF  is  fully  adaptive  because  all 
DOF  are  used  for  covariance  matrix  estimation.  Therefore,  twice  the  number  of  array 
elements  times  the  number  of  pulses  in  the  coherent  processing  interval  (CPI),  or 
2 NMP,  sample  support  vectors  are  required  to  operate  within  3  dB  of  optimal  when 
using  the  AMF. 
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2.3.3  FTS.  FTS  is  described  in  [5,13]  as  an  element-space  post-Doppler 
method  where  the  incoming  signal  from  each  antenna  element  is  Doppler  filtered  first. 
The  Doppler  filtered  data  is  then  used  to  estimate  the  interference  covariance  matrix 
for  spatially  adaptive  filtering.  A  window  is  typically  applied  across  Doppler  since 
FTS  has  no  temporal  degrees  of  freedom.  Post-Doppler  methods  such  as  FTS  intend 
to  localize  clutter  in  angle  in  order  to  reduce  the  adaptive  degrees  of  freedom.  FTS 
requires  2N  sample  support  vectors  to  obtain  performance  within  3  dB  of  optimal 
when  the  array  is  linear.  Two  separate  FTS  algorithms  are  presented  in  [5],  one  for 
uniformly  spaced  linear  arrays  and  one  for  uniformly  spaced  planar  arrays.  The  FTS 
algorithm  for  uniform  linear  arrays  is  shown  below  since  it  is  more  similar  to  the  FTS 
algorithm  derived  in  Chapter  III  for  non-uniformly  spaced  arrays.  The  Doppler  filter 
bank  is  defined  as 

F  =  diag(tb)W^  (2.25) 

where  W  is  an  M  x  M  discrete  Fourier  transform  (DFT)  matrix.  The  Mat  lab® 
function  “diag”  creates  a  diagonal  matrix  from  its  vector  argument  tb,  which  can  be 
a  window  function.  Otherwise,  tb  is  simply  a  vector  of  ones.  Equation  (2.25)  can  be 
written  as 

F=[f0  fi  ...  fM-i]  (2.26) 

showing  each  column  of  F  represents  a  different  Doppler  filter.  The  FTS  spatially 
adaptive  filter  weights  for  the  rnth  Doppler  bin  are 

W  m  =  Rm'a  (2.27) 

where 

Rm  =  (f m  <8>  Iv)i?R(fm  ®  Ijv)  (2.28) 

and  fm  is  the  mth  Doppler  filter  bin  of  F. 
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2.3.4  JDL.  The  JDL  algorithm  is  developed  in  [12]  and  presented  in  [5] 
for  the  uniformly  spaced  planar  array.  JDL  is  a  beamspace  post-Doppler  method 
where  adaptivity  is  applied  after  the  incoming  data  is  transformed  from  the  space- 
time  domain  to  the  anglc-Dopplcr  domain.  The  spatial  and  temporal  steering  vectors 
are  used  to  set  up  a  localized  processing  region  (LPR)  using  the  transformation  matrix 

T  =  [e(0t,  0_i)  e((j>t,  90)  0i)]  <g>  [b(d>_i)  b(w0)  b(u>i)] 

®[a(0_i ,0t)  a((j)0,6t)  a(0i,  0t)\  (2.29) 

where  90,  bo,  and  u)0  are  the  elevation  angle,  azimuth  angle,  and  Doppler  frequency 
at  the  LPR  center.  The  other  angles  and  Doppler  frequencies  are  simply  the  previ¬ 
ous  or  next  ones  in  the  LPR.  The  LPR  established  by  the  transformation  matrix  in 
Eqn.  (2.29)  has  3  degrees  of  freedom  in  azimuth,  3  degrees  of  freedom  in  Doppler, 
and  3  degrees  of  freedom  in  elevation.  The  azimuth,  Doppler,  and  elevation  degrees 
of  freedom  are  denoted  r]a,  rjb ,  and  ?/e,  respectively.  The  transformation  matrix  is  also 
used  to  transform  the  3D  steering  vector, 

v  =  Tffv.  (2.30) 

A  reduced  dimension  estimated  covariance  matrix  is  also  obtained  from  the  transfor¬ 
mation  matrix  as 

R  =  T^RT.  (2.31) 

The  filter  weights  are  then  computed  and  applied  based  on  the  transformed  3D  steer¬ 
ing  vector  and  estimated  covariance  matrix 

w  =  R  v 

=  (THRT)_1Tffv.  (2.32) 
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JDL  requires  2r]arji,r}e  sample  support  vectors  to  obtain  performance  within  3  dB  of 
optimal. 

2.4  Previous  Work  With  Circular  Arrays 

STAP  has  typically  been  used  with  uniformly  spaced  arrays  for  radar  systems 
in  the  past.  Non-uniformly  spaced  arrays  have  recently  been  studied  to  determine 
if  they  possess  any  advantages  over  uniform  arrays.  There  is  a  limited  amount  of 
information  regarding  STAP  performance  using  non-uniformly  spaced  arrays  in  radar 
systems.  The  articles  that  are  available  mostly  concentrate  on  circular  arrays,  partic¬ 
ularly  the  circular  array  researched  by  the  US  Navy’s  Ultra-High  Frequency  (UHF) 
Electronically  Scanned  Array  (UESA)  program.  A  circular  array  can  be  electronically 
scanned  across  360°  in  azimuth,  whereas  a  uniform  linear  array  must  be  mechanically 
steered  to  obtain  the  same  coverage. 

In  two  publications  by  Zatman  [14,15],  the  performance  of  STAP  for  airborne 
early  warning  (AEW)  radar  is  evaluated  when  applied  to  the  circular  UESA.  The 
circular  array  being  simulated  contains  60  elements.  However,  only  20  consecutive 
elements  are  active  when  transmitting  and  receiving.  Therefore,  the  antenna  can 
“rotate”  in  6°  increments  based  on  the  choice  of  excited  elements  and  scans  electron¬ 
ically  ±3°.  The  60  element  sidelooking  circular  array  is  compared  to  an  18  element 
sidelooking  uniform  linear  array.  The  two  arrays  have  the  same  overall  gain.  A  few 
other  specific  simulation  parameters  are  presented  in  Table  2.1.  However,  the  element 
spacing  is  not  given.  SINR  loss  is  used  as  a  performance  metric  to  compare  both  ar¬ 
rays.  Zatman’s  results  show  that  SINR  loss  does  not  change  as  a  function  of  range 
for  the  uniform  linear  array.  The  circular  array  has  a  wider  clutter  null,  a  higher 
clutter  rank,  and  has  about  0.5  dB  more  SINR  loss  than  the  uniform  linear  array. 
The  performance  of  each  array  is  also  evaluated  when  there  is  a  30°  array  rotation. 
Array  rotation  is  examined  to  study  effects  such  as  aircraft  crab  angle.  In  this  case, 
the  clutter  null  for  the  uniform  linear  array  is  broader  than  the  clutter  null  for  the 
circular  array.  Despite  the  wider  clutter  null,  the  uniform  linear  array  SINR  loss  is 
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Table  2.1:  Simulation  parameters  in  [14, 15]  for  18  element 

uniform  linear  array  and  60  element  circular  array. 


Variable 

Value 

M 

fo 

fr 

B 

ha  (aircraft  altitude) 
va  (aircraft  velocity) 
CNR 

Transmit  Taper 

18 

450  MHz 

300  Hz 

4  MHz 

9000  m 

100  knots 

45  dB  per  element  per  pulse  at  100  km 
Uniform  (None) 

still  about  0.5  dB  higher  than  the  circular  array  SINR  loss.  The  simulations  in  [14, 15] 
also  evaluate  pre-Doppler  and  post-Doppler  STAP  methods.  It  was  found  that  post- 
Doppler  STAP  performs  better  than  pre-Doppler  STAP  for  the  circular  array  since 
post-Doppler  STAP  reduces  clutter  rank  through  Doppler  filtering.  Range  varying 
adaptive  weights  can  be  used  to  mitigate  the  range  dependence  of  clutter  caused  by 
the  circular  array.  The  conclusions  by  Zatman  are  valuable  since  post-Doppler  STAP 
methods  are  applied  to  circular  arrays  in  Chapter  V. 

The  previously  established  mathematical  model  for  analyzing  an  airborne  radar 
system  with  a  sidelooking  uniform  array  discussed  in  this  chapter  provides  the  basis 
for  the  model  presented  in  the  following  chapter.  Chapter  III  changes  this  airborne 
radar  model  to  accommodate  non-uniformly  spaced  arrays.  A  radar  model  for  non- 
uniformly  spaced  arrays  is  valuable  because  any  array  configuration  can  be  analyzed. 
An  array  can  be  evaluated  to  see  how  detection  performance  is  affected  when  some  of 
its  elements  fail  to  operate.  Any  array  configuration  can  be  simulated  and  compared 
to  any  other  array  using  the  model  for  non-uniform  element  spacing.  This  model  is 
used  in  Chapter  V  to  simulate  circular  arrays  and  uniform  planar  arrays  with  element 
failure. 


16 


III.  Radar  Model  For  Non-Uniform  Element  Spacing 


This  chapter  incorporates  non-uniform  element  spacing  into  the  radar  data  model 
presented  in  [5].  Changing  the  element  positions  within  the  array  affects  the 
time/phase  delay  from  element-to-element.  This  new  time/phase  delay  causes  the 
spatial  steering  vectors  to  change,  bringing  about  changes  in  the  array  factor,  the 
space-time  snapshot  \i  and  covariance  matrices.  The  model  is  valid  for  side-looking 
linear  and  planar  arrays.  The  model  can  also  accommodate  3-dimensional  arrays. 
However,  shadowing  effects  have  not  been  taken  into  consideration.  Shadowing  occurs 
when  the  incoming  radar  wave  does  not  strike  some  array  elements  due  to  their 
location.  For  example,  this  may  occur  in  an  array  which  conforms  to  the  belly  of 
an  aircraft.  In  this  case,  the  elements  on  the  aircraft  underside  would  not  receive  a 
signal  coming  in  from  above  the  aircraft.  The  development  in  this  chapter  parallels 
the  development  presented  in  [5]  with  the  necessary  changes. 

3.1  Array  Geometry 

The  non-uniformly  spaced,  side-looking  array  being  modelled  now  contains  N 
elements  located  in  the  xz-plane  of  the  cartesian  coordinate  system.  The  radar  plat¬ 
form’s  velocity  vector  is  parallel  to  the  positive  rc-axis  shown  in  Fig.  3.1.  Any  element 
in  the  array  can  be  defined  as  the  reference  element.  The  reference  element  can  be 
located  anywhere  in  the  cartesian  coordinate  system  in  this  development.  An  incom¬ 
ing  radar  echo  strikes  the  reference  element  first  and  has  a  time  delay  associated  with 
all  other  elements.  The  reference  element  should  be  logically  chosen.  For  example, 
mathematical  evaluation  is  less  confusing  when  all  elements  arc  located  in  one  quad¬ 
rant  of  the  coordinate  system.  In  the  uniform  planar  array  situation,  it  makes  sense 
to  choose  the  upper  left  element  as  the  reference  element  as  in  Fig.  3.1.  The  array 
elements  can  also  be  numbered  in  any  order  and  the  model  will  still  be  valid. 

A  new  matrix  containing  the  arbitrary  element  positions  is  defined  as  D  which 
has  three  column  vectors  in  it  describing  the  x,  y,  and  z  positions  of  each  element 
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Figure  3.1:  Sidelooking  planar  array  lying  in  the  xz- plane. 

Boresight  is  defined  as  being  parallel  to  the  positive  y-axis.  The 
aircraft  velocity  vector  is  along  the  positive  x-axis. 

and  is  given  by 

D  =  [dx  dy  dz\.  (3.1) 

The  nth  row  of  D  describes  the  nth  element  location  and  is  written  as  D[n, :].  Using 
the  geometry  of  Fig.  3.1,  the  location  of  element  n  is 

D[n, :]  =  dx[rc]x  +  dy[n]y  +  dz[n]z,  (3.2) 

where  dx  [n] ,  dy[n],  and  dz[n]  are  scalar  distances  in  meters  measured  from  the  nth 
element  along  the  x-axis,  y-axis,  and  z-axis,  respectively,  and  x,  y,  and  z  are  unit 
vectors. 

It  is  common  to  describe  targets  relative  to  the  radar  using  a  coordinate  system 
other  than  the  cartesian  system.  The  radar  coordinate  system  shown  in  Fig.  3.2  is 
similar  to  the  spherical  coordinate  system.  Elevation  ( 9 )  and  azimuth  (0)  angles 
are  referenced  to  boresight  in  radar  coordinates.  Boresight  is  perpendicular  to  the 
array  face  and  is  along  the  positive  y-axis  in  this  development.  Positive  azimuth 
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Figure  3.2:  Sidelooking  planar  array  lying  in  the  xz- 
plane.  Negative  elevation  angle  is  measured  from  bore- 
sight  towards  the  ground  (negative  z-axis).  Positive  el¬ 
evation  angle  extends  from  boresight  towards  the  sky 
(positive  2- axis) .  Positive  azimuth  angle  is  measured 
from  boresight  towards  the  aircraft  nose  (positive  x- 
axis).  Negative  azimuth  angle  extends  from  boresight 
towards  the  aircraft  tail  (negative  re-axis). 

angles  are  measured  from  boresight  towards  the  aircraft  nose,  or  the  positive  re-axis, 
while  negative  azimuth  angles  are  measured  from  boresight  towards  the  tail.  Positive 
elevation  angles  are  measured  from  boresight  towards  the  sky,  or  the  positive  z-axis, 
while  negative  elevation  angles  are  measured  from  boresight  towards  the  ground.  It  is 
more  insightful  to  mathematically  analyze  the  radar  model  in  the  cartesian  coordinate 
system  since  it  is  the  most  familiar  coordinate  system.  Therefore,  a  transformation 
is  necessary  to  convert  an  azimuth  and  elevation  direction  in  the  radar  coordinate 
system  to  a  direction  in  cartesian  coordinates.  This  transformation  is  given  in  [13]  as 

k(0, 0)  =  cos  6  sin  0x  +  cos  6  cos  0y  +  sin  9 z.  (3.3) 
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A  new  element  position  matrix  D  is  now  set  np  to  accommodate  arbitrary 
array  configurations.  The  element  positions  affect  the  received  signal  time/phase 
delay  across  the  array.  The  transmit  signal  must  be  defined  before  the  received  signal 
can  be  described. 

3. 2  Transmit  Signal 

The  signal  transmitted  from  the  radar  is  mathematically  described  in  [13]  as 

s(t)  =  atu(t)e^Wot+lp\  (3.4) 

where  Eqn.  (3.4)  includes  a  random  phase  shift  tp,  a  sinusoidal  carrier  frequency  ca0, 
a  pulse  amplitude  at  (assumed  to  have  unity  energy),  and  an  envelope  function  u(t). 
The  pulse  width  and  pulse  repetition  interval  Tr  are  contained  in  the  envelope  function 
defined  in  [13]  as 

M— 1 

»«)  =  £  “p  (' t  —  mTr ).  (3.5) 

m= 0 

The  summation  of  M  pulses  defines  the  CPI,  which  is  a  block  of  data  available  for 
radar  processing.  The  radar  uses  this  block  of  data  to  statistically  model  the  environ¬ 
ment.  The  antenna  pattern  is  changed  based  on  the  CPI  to  null  out  interference  and 
determine  if  a  target  is  present.  The  received  signal  can  now  be  determined  since  it 
is  an  echo  response  of  the  transmitted  signal. 

3.3  Received  Signal 

The  signal  received  at  the  nth  array  element  has  the  same  form  as  the  transmit¬ 
ted  signal  except  that  it  is  delayed  in  time  by  r  and  shifted  in  frequency  by  ft  due  to 
the  target  distance  and  velocity,  respectively.  The  signal  delay  equals  the  round  trip 
time  to  the  target  Tri  plus  the  time  it  takes  the  signal  to  travel  from  the  reference 
element  across  the  array  to  the  nth  element  r'n.  The  complete  signal  delay  is 

Tn  =  Tr  +  t' .  (3.6) 
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The  signal  received  at  the  nth  element  is  then  written  as 


s 


n 


s(t  -  Tn)ej2irMt-Tn) 

aru(t  -  Tn)ej 2*Mt-Tn)eV*Mt-Tn)ej'P' 


(3.7) 


The  received  pulse  amplitude  ar  is  different  than  the  transmitted  pulse  amplitude  at 
due  to  various  factors  such  as  the  target’s  radar  cross  section  (RCS)  a,  atmospheric 
attenuation,  system  losses  LS)  etc.  as  encapsulated  in  the  radar  equation  of  [10], 

_  PtGAea 

max  ~  (4tt )*kT0BFn{S/N)MLs 


(3.8) 


Pt  is  transmit  power,  G  is  transmit  antenna  gain,  Ae  is  the  effective  received  antenna 
aperture,  k  is  Boltzmann’s  constant  1.38  x  10-23  J/deg,  T0  is  standard  temperature 
290  K,  B  is  receiver  bandwidth,  Fn  is  receiver  noise  figure,  and  ( S/N ) M  is  the  effective 
SNR  of  M  integrated  pulses. 

The  time  delay  from  the  reference  element  to  the  nth  element  is  found  by  pro¬ 
jecting  the  target  direction  unit  vector  onto  a  vector  representing  the  nth  element’s 
location, 

/  _  k(0,0)»D[n,:] 

Tn  c 

(cos  9  sin  0x  +  cos  6  cos  4>y  +  sin  9z)  •  (d-Jnjx  +  d.y[n]y  +  d2[n]z) 

c 

d,j.  [n]  cos  9  sin  0  +  dy  [n]  cos  9  cos  (j)  +  dz  [n]  sin  9 


Then,  the  phase  delay  to  the  nth  element  is  simply  the  corresponding  time  delay 
multiplied  by  c uQ, 


uaTn  =  2vr/0 


d^.  [n]  cos  9  sin  cj)  +  dy  [n]  cos  9  cos  (p  +  dz  [n]  sin  9 


(3.10) 
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When  Eqn.  (3.6)  is  substituted  into  Eqn.  (3.7),  the  received  signal  becomes 

sn(t)  =  aru(t  -Tr-  T^ej2n{fo+m-TR-T^e^.  (3.16) 

A  few  assumptions  can  be  made  to  reduce  Eqn.  (3.16).  The  first  is  the  narrowband 
assumption:  r'n  <C  pulse  width.  In  other  words,  the  time  delay  from  the  reference 
element  to  the  nth  element  can  be  neglected  in  the  envelope  function  since  it  is  much 
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smaller  than  the  pulse  width.  Equation  (3.16)  now  becomes 


sn(t)  =  aru(t  -  TR)eJ’27r(/o+/t)te-j27r(/o+/t)Tile-:''2T(/o+/t)r"e^.  (3.17) 

The  spatial  frequencies  of  Eqns.  (3.12)  through  (3.14)  can  be  substituted  into  the 
third  exponential  term  of  the  received  signal 

e~j2irf0Tll  _  e-ju0T'n  _  e-j2n(&x[n]+'dy[n}+'dz[n})  (3.18) 


yielding 

sn{t)  =  aru(t  -  TR)e^(/o+/t)te-j2^{/o+/t)Tfle^2^/tT™e^2T(^[n]+^w+,9z[n])e^.  (3.19) 

Examining  the  received  signal  again  leads  to  another  reduction.  The  second  expo¬ 
nential  term  is  constant,  so  it  can  be  grouped  into  the  arbitrary  phase  term  e^.  The 
third  exponential  term  can  be  approximated  as  1  because  the  target  Doppler  shift  ft 
is  assumed  to  be  small  relative  to  the  carrier  frequency  fQ  and  the  time  delay  r'n  has 
the  speed-of-light  in  its  denominator.  The  received  signal  is  now  represented  as 

sn(t)  =  aru(t  -  TR)ej2^(/o+/t)te-i2T(^[n]+^[nl+^[nl)ejV.  (3.20) 

After  the  signal  is  received,  it  is  frequency  down-converted  by  a  mixer  for  pro¬ 
cessing.  Figure  3.3  shows  a  simple  block  diagram  of  the  receiver  chain  for  one  an¬ 
tenna  element.  The  mixing  process  converts  the  signal  frequency  to  an  intermediate 
frequency  (IF).  In  this  derivation,  the  IF  is  assumed  to  be  the  baseband  frequency. 
Therefore,  the  carrier  frequency  is  no  longer  in  the  mixed  signal 

sn(t)  =  aru(t  -  riJ)ei2^tte-j2"(^N+^[nl+^[n])e^.  (3.21) 

After  mixing,  the  received  signal  is  match  filtered  to  maximize  echoes  from  the  desired 
radar  signal  and  to  attenuate  unwanted  signals.  The  matched  filter  output  is  obtained 
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Antenna 

Element 


xH(k) 


Figure  3.3:  Basic  block  diagram  of  receiver  for  one  antenna 
element.  The  signal  is  received  by  the  antenna  element.  Then, 
the  frequency  is  down  converted  to  an  intermediate  frequency 
by  a  mixer.  After  mixing,  the  signal  is  match  filtered.  The 
signal  is  then  sampled  by  an  analog-to-digital  converter  for 
processing. 

by  convolving  the  filter  impulse  response  with  the  received  signal, 


xn{t) 


sn{S)h(t  —  5)d8. 


(3.22) 


The  filter  impulse  response  is 

h(t)  =  u*p{-t )  (3.23) 

since  the  received  signal  is  match  filtered  pulse-by-pulse.  After  substituting  the  re¬ 
ceived  signal  and  filter  impulse  response  into  Eqn.  (3.22),  the  matched  filter  output 
becomes 


xn{t) 


a  M  +  &y  [n] [ra] )  ejip 


M—l 

E  up -  Tr  - 

m= 0 


mTr)u;(8  -  t)ej2nftSd8. 


(3.24) 


(3.25) 


The  matched  filter  output  equation  can  be  reduced  by  using  a  change  of  variables 
within  the  integral.  Let  /3  =  8  —  Tr  —  mTr  =>■  8  —  (3  +  Tr  +  mTr  and  df3/d8  =  1.  Now 
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the  matched  filter  output  is 


M—l 


Xn(t)  =  are-32^(^x[n]+'&yln]+‘&z[n])ejipej2wftTR  ej2-KftmTr 


m= 0 


up((3)u*J(3  +  Tr  +  mTr  -  t)ej2nMd(5 .  (3.26) 


The  integral  term  in  Eqn.  (3.26),  known  as  the  time  frequency  autocorrelation  function 
(TFACF)  for  a  transmitted  pulse,  can  be  approximated  as  1  because  the  transmitted 
pulse  was  assumed  to  have  unity  energy  and  the  rnth  pulse  has  a  range  delay  of 
Tr  +  mTr  at  the  range  cell  under  test.  This  translates  into  a  perfect  match  in  range 
at  the  TFACF  peak.  The  waveform  is  also  assumed  to  be  somewhat  Doppler  tolerant. 
Therefore,  the  match  filtered  signal  for  the  array’s  nth  element  and  for  a  target  at  some 
range  is 

M—l 

_  areJ‘Pe-j2n(&x[n]+'&y[n]+&z[n])  E  2tt nuu  27) 

m= 0 

where  u  =  ft/fr  is  the  normalized  Doppler  frequency  and  fr  is  the  pulse  repetition 
frequency 

fr  =  1  /Tr.  (3.28) 


If  the  normalized  Doppler  frequency  is  greater  than  one,  then  the  corresponding  ve¬ 
locity  is  ambiguous  to  the  radar.  The  e:’2n^tTR  term  is  lumped  into  the  random  phase 
term  since  it  is  constant.  The  mth  match  filtered  pulse  for  the  nth  element  can  be 
expressed  as 


—  m  p-j2ir(^*H+Oj/[n]+tfz[n])_;27rma; 
nm  °  i 


(3.29) 


where  at  is  a  complex  value  representing  the  amplitude  and  random  phase  of  the 
pulse.  The  match  filtered  signal  is  now  in  a  form  representing  the  response  at  each 
antenna  element  due  to  each  pulse.  From  here,  it  is  necessary  to  develop  an  expression 
representing  the  returns  from  all  pulses  across  the  entire  array  because  the  radar  uses 
all  the  data  within  a  CPI  for  target  detection.  First,  the  returns  due  to  the  nth  element 
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and  all  M  pulses  can  be  written  in  an  M  x  1  vector  as 


;n  =  ate-j2^[n]e-j2n®y[n]e-j2n^[n]h(0) 


(3.30) 


where 

b(u>)  =  [1  ej2™  ...  epMM-^T  (3.31) 

is  the  temporal  steering  vector  dehned  in  [13].  Equation  (3.30)  is  now  a  relaxed  form 
of  the  incoming  returns  dehned  in  [5]  and  accommodates  arbitrary  element  spacing. 
At  this  point,  there  are  a  few  different  approaches  that  can  be  taken  to  represent  the 
phase  progression  across  all  elements.  One  way  is  to  write  a  column  vector  containing 
the  phase  progression  to  each  element.  Writing  such  a  vector  is  the  most  direct 
method  although  it  may  be  tedious  for  arrays  with  many  elements.  This  vector  of 
element  phase  progressions  is  called  c  and  is  written  as 

C(Q  0)  _  [e-j2n0x[O}e-j2w&y[O]e-j2ir‘&z[O]  e-j2n‘&x[i\e-j2ir0y[i\e-j2ir0z[l] 

g  —j2ir&x [JV— 1] g —j2n‘0y [JV—  1] g  —j2ir&z [TV— 1] j T  (3  32) 

The  Kronecker  Product  (see  Appendix  A)  can  then  be  taken  with  the  c  vector  and 
the  temporal  steering  vector  to  produce  a  vector  known  as  the  space-time  snapshot 
X •  The  space-time  snapshot  contains  the  returns  from  all  elements  and  pulses  and  is 


Xt  =  «(bH®c  (0,0) 


(3.33) 


for  a  single  point  scatterer. 
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Expanding  Xt  results  in 


e-j2ir&x  [0]  e~j27rdy  [0]  g-j27ri9z  [0] 
e~j2ir&x  [1]  e-j2irOy  [1]  e~j2ir&z  [1] 

e— j2ir&x  [N- 1]  g  —  j 2ir&y  [N— 1]  g  -  j 2tti?z  [AT— 1] 
gj2irtD  g  -  j2ir&x  [0]  g  ■ -  j2ir&y  [0]  g  -  j  2tti?z  [0] 
ej2irui e~j2TT'3x  [l]e-j2ir&y  [l]e-j2w&z  [1] 

Xt  =  at  |  : 

gj27ra>  e-j2ir&x  [JV — 1]  e-j2ir0y  [JV— 1]  e~j27rdz  [N— 

e  j27r  (M-  l)a)  g- j27ri9x  [0]  g  -  j  2^1?^  [0]  g  -  j  2?n9z  [0] 
gj'27r(M-l)a)g-j27ri9x[l]g-j27ri?!/[l]g-j27ri?z[l] 

gj27r(M-l)a)g-j27ri>x  [TV-l]g-j27ri9H[Af-l]g-j27ri9z[ 

The  space-time  snapshot  contains  the  target  amplitude  and  direction.  The  target 
spatial  location  is  represented  by  the  3D  steering  vector 

v(u>,  6,  4>)  =  b(o>)  (8)  c(9, 4>).  (3.34) 

The  3D  steering  vector  describes  the  specific  azimuth  angle,  elevation  angle,  and 
normalized  Doppler  frequency  of  a  target  at  a  particular  location.  The  target  return 
dependence  on  the  3D  steering  vector  can  now  be  seen  when  Eqn.  (3.34)  is  substituted 
into  Eqn.  (3.33),  resulting  in 

Xt  =  atv(Q,9,(f>).  (3.35) 

In  other  words,  the  space-time  snapshot  is  simply  the  3D  steering  vector  scaled  in 
amplitude  and  phase. 

The  3D  steering  vector  can  be  used  to  maximize  the  radar  response  in  the  radar 
“look  direction.”  The  complex  conjugate  of  the  3D  steering  vector  must  be  taken  to 
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maximize  the  radar  response.  This  complex  conjugated  steering  vector  can  be  thought 
of  as  the  steering  vector  from  the  target  to  the  radar,  representing  the  expected  phase 
delays  across  the  array  from  a  particular  location  in  space.  These  expected  phase 
delays  are  essentially  filter  weights  for  each  element  and  pulse  allowing  target  returns 
to  be  constructively  added.  The  filter  weight  vector  is 

w  =  v(o>,  6,  0).  (3.36) 

This  particular  weight  vector  is  non-adaptive  because  it  simply  represents  the  ex¬ 
pected  phase  delays  from  a  target  at  a  specific  location  in  space  without  considering 
the  interference  from  the  surrounding  environment.  The  filter  weights  are  mathemat¬ 
ically  applied  to  the  target  returns  by  the  inner  product  between  the  filter  weight 
vector  and  the  space-time  snapshot,  producing  an  output  given  by 

Output  =  wHxt  = 

=  o;iV^(d),6>»v(d),6>,0) 

N—l  M—l 

=  «<££  ej2n('ax  [n]+^y[n]+‘0z[n\)e—j2'Kmuie—j2'K('&x[n\+'Oy[n]+‘dz[n])ej2'Kmu> 

n= 0  m= 0 

=  atMN.  (3.37) 

The  multiplication  and  addition  between  the  filter  weights  and  the  space-time  snap¬ 
shot  represent  coherent  integration  over  M  pulses  across  the  entire  array.  The  N 
factor  in  Eqn.  (3.37)  is  the  spatial  coherent  integration  gain  across  the  array,  while 
M  is  the  coherent  integration  gain  from  integrating  M  pulses. 

The  spatial  steering  vector  is  now  defined  for  arbitrary  element  spacing,  allowing 
the  3D  steering  vector  and  space-time  snapshot  for  a  scatterer  to  be  determined.  The 
spatial  steering  vector  is  used  to  determine  the  array  pattern,  jammer  model,  and 
clutter  model.  The  array  pattern  is  useful  because  it  shows  how  the  array  illuminates 
the  environment.  The  element  pattern  must  be  defined  before  the  array  pattern  is 
discussed. 


3.4  Element  Pattern 


Each  element  in  the  array  is  assumed  to  have  identical  voltage  and  power  pat¬ 
terns  as  discussed  in  [6, 13].  The  voltage  pattern  as  defined  for  a  planar  array  in  [5] 

is 

f  008(6*)  cos(0)  —90°  <  (f>,9  <  90° 

/(M)  =  \ 

be  cos (9)  cos(</>)  90°  >  0,  9  <  270° 

where  be  is  a  backlobe  scaling  factor.  Mutual  coupling  effects  are  not  considered  in 
this  element  pattern  expression.  The  voltage  pattern  is  described  as  a  cosine  pattern 
with  backlobe  attenuation  and  is  shown  in  Fig.  3.4.  The  backlobe  factor  be  models  the 
attenuation  resulting  from  the  surfaces  behind  the  antenna  elements.  Such  surfaces 
may  include  the  back  of  the  array,  the  aircraft  fuselage,  or  anything  else  behind  the 
array  that  attenuates  the  transmitted  signal.  The  element  power  pattern  equals  the 


Figure  3.4:  3D  element  voltage  pattern  converted 

from  radar  coordinates  to  cartesian  coordinates.  This 
plot  illustrates  the  energy  emitted  from  each  antenna 
element  in  a  given  direction.  The  antenna  element  is 
located  at  the  origin  (x  —  0,  y  —  0,  z  —  0)  for  this 
illustration. 
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magnitude  of  the  element  voltage  pattern  squared 


g{QA)  =  l/(M)|2-  (3-38) 

The  array  pattern  is  defined  in  the  next  section  since  it  is  based  on  the  element 
pattern. 

3.5  Array  Pattern 

Once  the  voltage  pattern  is  established  for  each  element,  the  array  pattern  must 
be  determined.  The  array  pattern  represents  the  overall  transmitted  voltage  pattern 
from  the  entire  array.  Sidelobes,  nulls,  and  backlobes  are  illustrated  in  the  array 
pattern.  The  array  pattern  is  constructed  from  the  array  factor  and  the  element 
pattern.  The  array  factor  is  the  sum  of  all  the  phase  weights  across  the  array  and  is 
expressed  as 

w(e,cf>)  =  c(e,  ct>)Hc(e',  t'),  (3.39) 

where  vector  c (9,  <fi)  contains  the  target  location  in  space  and  vector  c (O',  <fi')  contains 
the  desired  mainbeam  location.  Array  beamforming  theory  is  discussed  more  thor¬ 
oughly  in  [1,11].  The  array  factor  is  simply  scaled  by  the  element  pattern  to  obtain 
the  array  pattern.  The  array  voltage  pattern  is 

F(M)  =  inM)/(M)  (3.40) 

and  the  array  power  pattern  is 

G(M)  =  \w{9,<t>)f{0,<t>)\2  =  |WW)I2$(M).  (3.41) 

A  steerable  antenna  pattern  can  now  be  generated  for  an  array  with  arbitrary 
element  locations.  The  noise,  jammer,  and  clutter  models  are  described  next  because 
they  are  directly  affected  by  the  antenna  pattern. 
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3. 6  Interference  Models 

The  next  step  in  updating  the  radar  data  model  is  to  reevaluate  the  interference 
models.  Covariance  matrices  are  used  in  the  interference  models  to  describe  the 
statistical  expectations  of  power  across  all  elements  and  pulses  at  a  specific  range  cell 
for  a  given  source  of  interference.  It  is  important  to  note  that  the  noise,  jammer,  and 
clutter  models  are  not  actually  being  changed  here.  These  models  are  simply  being 
modified  so  interference  can  be  modelled  for  non-uniformly  spaced  arrays. 

3. 6. 1  Noise  Model.  Within  the  radar  system,  each  array  element  has  its  own 
receiver.  As  a  result,  noise  power  is  not  correlated  from  element-to-element.  Noise 
power  is  also  uncorrelated  from  pulse-to-pulse  since  it  is  assumed  to  be  white  noise. 
This  assumption  is  only  applicable  when  the  PRF  is  much  less  than  the  waveform 
bandwidth  and  the  waveform  bandwidth  is  much  less  than  the  carrier  frequency.  The 
resulting  expected  noise  power  for  any  given  element  and  pulse  is  very  similar  to 
the  expected  value  of  any  other  element  and  pulse.  Therefore,  the  noise  covariance 
matrix  established  in  [13]  is  the  noise  power  per  element  per  pulse  o'2  multiplied  by 
an  identity  matrix  of  size  MN 

R«  =  e{xnXn  }  =  (3.42) 


!  The  noise  covariance  matrix  in  [13]  is  for  a  linear  array  containing  N 
elements.  The  noise  covariance  matrix  in  Eqn.  (3.42)  is  for  a  linear,  planar, 
or  3D  array  with  N  total  elements. 

It  can  be  seen  from  Eqn.  (3.42)  that  the  noise  covariance  matrix  only  depends  on 
the  number  of  elements  and  pulses.  Element  position  on  the  array  does  not  affect  the 
expected  noise  power.  Therefore,  the  noise  covariance  matrix  is  identical  for  arrays 
with  the  same  number  of  elements  despite  the  element  positions. 

3.6.2  Barrage  Noise  Jammer  Model.  The  jammer  model  is  for  a  barrage 
noise  jammer  [5,6,13].  Although  it  is  not  used  in  simulations  in  the  following  chapters, 
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the  jammer  model  is  presented  here  for  completeness.  Its  power  spectral  density 


received  at  any  element  is 

%(M)  Ag 

°  (4vr)2i?fLr  ‘ 

Lr  is  receiver  losses,  Rj  is  the  distance  between  the  radar  and  jammer,  and  S3  is  the 
jammer’s  effective  transmitted  power  spectral  density.  The  expected  jammer  power 
is  spatially  correlated  since  the  jammer  is  located  at  a  specific  elevation  angle  and 
azimuth  angle.  However,  there  is  no  temporal  correlation  because  the  transmitted 
jammer  amplitude  is  complex  and  varies  from  pulse-to-pulse.  The  random  jammer 
amplitudes  are  written  in  vector  form  as 

ctj  =  [tt0  aL  ...  aM- i]T-  (3.44) 

The  size  of  this  vector  is  M  x  1  since  there  is  a  different  amplitude  for  each  pulse. 
The  space-time  snapshot  for  the  barrage  noise  jammer  is  expressed  as 


Xj  =  aj®  c(0,0).  (3.45) 

The  jammer  covariance  matrix  can  be  found  from  the  space-time  snapshot  as 

Rj  =  £{XjXf}  =  £{[<*]  ®  c(0,0)][aj  ®c(d:(p)}H}.  (3.46) 

If  the  transmitted  jammer  power  is  assumed  to  be  stationary  over  a  CPI,  then 

£{a.jCtf}  =  0-%-Im,  (3.47) 

where 

6  -  Jt  <3-«> 

is  the  jammer-to-noise  ratio  per  element.  This  assumption  reduces  the  jammer  co- 
variance  matrix  to 

Rj  =  a2ijlM  ®  c(0,  0)c(0,  (f))H .  (3.49) 
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3.6.3  Clutter  Model.  The  clutter  covariance  matrix  describes  the  natural 
environment  surrounding  the  platform.  Radar  signals  bounce  off  trees,  buildings, 
and  other  objects  in  the  environment.  The  clutter  model  used  is  a  constant  gamma 
model  [10]  for  ground  or  sea  clutter.  Sky  clutter  is  neglected  in  this  model  because  it 
is  so  small  compared  to  ground  clutter.  A  few  different  distances  and  angles  defined 
in  [5]  are  reiterated  here  for  the  clutter  model  development.  The  earth  is  considered 
to  be  a  perfect  sphere  with  an  effective  radius  4/3  its  normal  size,  ae  =  4re/3.  Clutter 
is  divided  into  rings  and  patches  around  the  aircraft  as  in  Fig.  3.5.  Each  clutter  ring 
surrounds  the  aircraft  at  a  range  Rc.  The  elevation  angle  from  the  aircraft  to  the 
clutter  ring  is 


6r  =  —  sin 


-l 


v  R2c  +  ha(ha  +  2  ae) 


(3.50) 


2  Rc(ae  +  ha ) 

where  ha  is  the  aircraft  altitude.  Each  clutter  ring  is  divided  into  Nc  patches.  The 


Figure  3.5:  Overhead  view  of  clutter  rings  and  patches 
surrounding  the  radar  platform. 
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grazing  angle  i^c  shown  in  Fig.  3.6  is  defined  as  the  angle  between  a  line  tangent  to 
the  earth  at  the  clutter  patch  and  a  line  from  the  radar  to  the  clutter  patch.  The 
grazing  angle  is  determined  as 


-0c  =  —  sin 


ha(ha  T  2ae) 

‘2Rco,e 


(3.51) 


The  grazing  angle  is  set  equal  to  0  to  determine  the  range  from  the  aircraft  to  the 
horizon.  The  horizon  range  is 


Rh  \J h ^  T  2haae. 


(3.52) 


The  unambiguous  range  Ru  =  cTr/ 2  is  used  with  the  horizon  range  to  determine  the 
number  of  range  ambiguities  Nr  as 


v 


RcNr  <  Rh 

RcNr  >  Rh, 


where  ("x]  means  round  up  to  the  nearest  integer  and  |ycj  means  round  down  to 


Figure  3.6:  Side  view  of  clutter  ring  layout,  showing  the  grazing  angle. 
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the  nearest  integer.  The  ambiguous  range  rings  are  at  Ri  =  iRu  +  Rc  where  i  = 

0, 1,  ...,Nr  —  1. 

The  major  change  that  takes  place  in  the  clutter  model  for  non-uniform  element 
spacing  is  the  spatial  frequency  of  the  ikth  clutter  patch,  where  the  i  subscript  defines 
the  current  clutter  ring  based  on  elevation  angle  6i  and  the  k  subscript  designates 
the  current  patch  around  the  clutter  ring  based  on  azimuth  angle  (f)k-  The  vector 
describing  the  elevation  and  azimuth  angle  of  the  ikth  clutter  patch  is  projected  onto 
the  array’s  nth  element  and  divided  by  the  wavelength.  The  spatial  frequency  defined 
in  Eqn.  (3.11)  is  evaluated  at  0*  and  (f)k 


&ik[n. 


k (9i,<t>k)  •  D[n, : 

A0 

dx  [n]  cos  6 1  sin  <i>k  +  dy  [ n ]  cos  0*  cos  (t>k  +  d2  [n]  sin  6i 


(3.53) 


for  the  ikth  clutter  patch.  The  Doppler  frequency  of  the  ikth  clutter  patch  must  also 
be  determined  in  order  to  construct  the  clutter  space-time  snapshot.  Using  Eqn.  (3.3), 
the  ikth  clutter  patch  Doppler  frequency  is  defined  as 


fc(0i,<t>k ) 


2k (0j,  (pk)  •  va 

A0 

2va  cos  6j  sin  (f>k 


(3.54) 


where  va  =  uax  is  the  aircraft  velocity  vector.  The  clutter  space-time  snapshot  is 
determined  as 

Nr  —  1  Nc  —  1 


xc=  aikb(ujik)  (8)  c (9i,<j>k), 


(3.55) 


i= 0  fc=0 

using  the  normalized  Doppler  frequency  Olk,  elevation  angle  6t,  and  azimuth  angle  (f)k 
of  the  ikth  clutter  patch.  The  summations  in  Eqn.  (3.55)  physically  represent  the  sum 
of  all  returns  from  each  clutter  patch  around  each  range  ring.  The  random  clutter 
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patch  amplitude  is  found  based  on  its  RCS.  The  ikth  clutter  patch’s  RCS  is 


Cifc  =  cro(6>j,0jfc)i?iA0ARsec(V'j),  (3.56) 

where  A cj)  =  ( 2ir/Nc )  and  A R  =  c/(2B).  The  radar  system  bandwidth  is  B.  Each 
clutter  patch  is  assumed  to  have  a  specific  backscatter  reflectivity  described  by 


cro(0ir<f>k)  =  7  sbi^i, 


(3.57) 


where  7/y  is  the  grazing  angle.  The  7  term  characterizes  clutter.  A  worst  case  clutter 
environment  may  have  a  7  value  as  high  as  -3  dB,  suggesting  the  clutter  is  very  strong 
and  may  be  caused  by  an  urban  environment  [8].  Clutter  can  be  as  weak  as  -29  dB, 
indicating  flat  terrain  such  as  a  desert  environment  [8].  After  the  clutter  patch  RCS 
is  determined,  the  element  clutter-to-noise  ratio  (CNR)  can  be  calculated  as 


Cik 


PtGt(0j,  4>k)g{0j,  <t>k)Xjaik 

{4n)3N0BLsRf 


(3.58) 


using  the  radar  equation.  The  clutter  patch  amplitude  is  found  by  simply  multiplying 
the  CNR  by  the  noise  power  a'2. 

Once  the  clutter  patch  amplitude  is  found,  the  clutter  space-time  snapshot  can 
be  computed  using  Eqn.  (3.55).  Then,  the  clutter  covariance  matrix  is  found  by  taking 
the  expected  value  of  the  space-time  snapshot’s  outer  product 


Nr-l  Nc- 1 

R-c  =  £{XcXc  }  =  0-2  &kHuik)Huik)H  ®  C (0i,  </>fc)c(^>  (pk)H ■  (3.59) 

i= 0  k= 0 


The  radar  model  is  now  in  a  format  suitable  for  radar  system  modelling  and 
analysis  of  non-uniformly  spaced  arrays.  Antenna  patterns  can  be  generated  and 
interference  can  be  modelled.  The  next  logical  step  is  to  update  the  STAP  models  so 
STAP  can  be  applied  to  non-uniform  arrays. 
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3.7  STAP  Models 


Two  practical  STAP  algorithms,  JDL  and  FTS,  are  updated  in  this  section. 
The  MF  and  AMF  algorithms  discussed  in  Chapter  II  do  not  need  to  be  modified 
because  the  3D  steering  vector  v  accounts  for  the  non-uniform  spacing. 

3. 7. 1  Joint  Domain  Localized.  The  JDL  model  changes  based  on  the  spatial 
steering  vector  c (6,  4>).  The  incoming  space-time  snapshot  data  is  transformed  to  the 
angle-Dopplcr  domain  using  the  redefined  transformation  matrix 

T  =  [b(d>_i)  b(cJo)  b(wi)]  <g>  [c(0_i,0_i)  c(0_i,0o)  c(0_ i,0i) 

c(0o,0_i)  c(0o,  00 )  c(0o,0i)  c(0!,0_i)  c(0i,0o)  c(0!,0i)].  (3.60) 

The  spatial  steering  vector  c (0,  (j) )  in  Eqn.  (3.60)  contains  the  phase  progression  to 
each  element  based  on  the  azimuth  and  elevation  angle  of  interest,  4>o  and  6q.  The 
resulting  transformation  matrix  is  of  size  NM  x  rp,'qa;qe  where  r/ft,  qa,  and  qe  are  the 
degrees  of  freedom  in  Doppler,  azimuth,  and  elevation,  respectively.  T  in  Eqn.  (3.60) 
has  3x3x3  degrees  of  freedom  since  there  are  3  different  Doppler  values,  3  different 
azimuth  angles,  and  3  different  elevation  angles.  This  3x3x3  region  contained  within 
the  entire  angle-Dopplcr  space  is  known  as  the  localized  processing  region  (LPR). 
Once  the  transformation  matrix  has  been  updated,  the  rest  of  the  JDL  algorithm  is 
identical  to  the  one  presented  in  [5]  and  discussed  in  Chapter  II.  The  transformed 
space-time  snapshot  and  steering  vector  are  still  x  —  and  v  =  T^v,  respectively. 
The  estimated  covariance  matrix  is  R  =  e{xXH }  =  T^RT.  The  adaptive  filter  is 
w  =  R  V. 

3.7.2  Factored  Time- Space.  The  FTS  algorithm  in  [5,13]  also  has  a  slight 
modification  based  on  the  new  radar  model  for  non-uniform  element  spacing.  The 
Doppler  filter  bank  is  the  same  as  in  [5], 

F  =  diag(tb)W^.  (3.61) 
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The  “diag”  term  is  a  Mat  lab®  function  that  creates  a  diagonal  matrix  from  the  vector 
argument  and  XVH  is  an  M  x  M  DFT  matrix.  The  tb  term  can  either  be  a  vector  of 
ones  or  a  window  function.  If  tb  is  a  vector  of  ones,  then  each  column  of  F  is  simply 
the  temporal  steering  vector  b(u;)  for  a  different  Doppler  frequency.  Equation  (3.61) 
can  be  expanded  as 

F=[f0  fi  ...  (3.62) 

showing  each  column  as  the  Doppler  filter  for  each  bin.  The  output  from  Doppler  bin 
m  is 

Xm  =  (frn  ®  In)H  X  (3.63) 

when  applied  to  the  space-time  snapshot.  The  spatially  adaptive  filter  weights  for  the 
mth  Doppler  bin  are 

wm  =  Rm1c(0,0),  (3.64) 

where  the  covariance  matrix 

R-m  =  s{xmXm}  =  (fm  ®  Ijv)^  R  (f, m  ®  I n)  (3-65) 

is  now  of  size  N  x  N,  reducing  the  required  sample  support.  The  rest  of  the  FTS 
algorithm  is  the  same  as  [5].  For  example,  the  full  dimensioned  weight  vector  is 

wm  =  fm  0  wm.  (3.66) 


3. 8  Summary 

The  radar  model  has  been  updated  to  incorporate  non-uniform  element  spacing 
for  linear  and  planar  arrays.  The  changes  began  by  redefining  the  position  vector 
describing  each  element  location  in  the  array.  The  new  element  position  definition 
affects  the  propagating  wave  time/phase  delay  across  the  array,  the  spatial  frequencies, 
the  steering  vectors,  the  array  pattern,  the  space-time  snapshots,  and  the  covariance 
matrices.  The  JDL  and  FTS  STAP  algorithms  have  also  been  modified  so  they  are 
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compatible  with  the  new  airborne  radar  model  for  non-uniform  element  spacing.  The 
model  is  in  a  format  suitable  for  data  generation  and  analysis. 
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IV.  Radar  Model  Verification 


The  new  radar  model  derived  in  Chapter  III  for  non-uniform  element  spacing 
is  verified  in  this  chapter.  Verification  is  accomplished  by  inserting  uniform 
element  spacing  into  the  model  derived  for  non-uniform  element  spacing.  Then,  the 
radar  model  presented  in  Chapter  III  is  compared  to  the  models  in  [13]  and  [5]  to  see 
if  they  produce  the  same  results.  The  model  presented  in  Chapter  III  will  identically 
collapse  back  to  the  model  in  [5, 13]  only  when  the  array  is  linear.  The  two  models 
still  produce  the  same  results  when  the  array  is  planar,  except  the  steering  vectors, 
the  space-time  snapshots,  and  the  covariance  matrices  are  ordered  differently  in  each 
model.  Therefore,  Matlab®  simulation  is  used  to  verify  that  the  model  in  Chapter  III 
produces  the  same  results  as  the  model  in  [5]  when  the  planar  array  contains  uniform 
element  spacing. 

4-1  Linear  Array 

The  model  presented  in  [5]  is  for  a  planar  array  with  N  azimuth  channels  and 
P  elevation  channels.  When  P  =  1,  the  array  contains  N  elements  extending  along 
the  negative  re-axis  with  the  reference  element  located  at  the  origin,  as  in  [13].  The 
inter-element  spacing  between  each  element  is  dx.  The  position  vector  of  each  element 
defined  in  Chapter  III  as 

D[0, :]  =  dx[0]x  +  dy[0]y  +  d2[0]z 
D[l, :]  =  dx[l]x  +  d.y[l]y  +  d2[l]z 

D [N  -  2, :]  =  d X[N  -  2]x  +  d y[N  -  2]y  +  d2[V  -  2]z 
D [N  -  1, :]  =  d X[N  -  l]x  +  d y[N  -  l]y  +  d2[lV  -  l]z 
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can  now  be  written  as 


D[0, :]  =  -(0 )dxx  +  (0)y  +  (0)z 
D[l, :]  =  -(1  )dxx  +  (0)y  +  (0)z 


D [N  -  2, :]  =  —  (JV  -  2 )dxx  +  (0)y  +  (0)z 
D[N  -  1, :]  =  -(N  -  1  )4x  +  (0)y  +  (0)z 

due  to  uniform  element  spacing.  By  examining  the  above  element  position  vectors, 
the  nth  element  location  can  be  expressed  as 

D[n, :]  =  —  (n)dxit  +  Oy  +  Oz 

=  — ndTx.  (4.1) 


The  y  and  z  terms  are  neglected  since  the  linear  array  is  located  along  the  x-axis  at 
dy  —  dz  =  0.  The  position  vector  is  consistent  with  [5]  for  a  linear  array  and  flows 
through  the  entire  radar  model.  For  example,  the  spatial  frequency  becomes 


i9[n] 


k(6*,  (p )  •  D[n, : 
dx  [n]  cos  6  sin  0 


— ndx  cos  9  sin  (p 


(4.2) 


The  spatial  frequency  ^-component  is  defined  in  [13]  as  a  scalar 


9t  = 


dx  cos  9  sin  0 


(4.3) 
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so  that  1 9[n]  =  —  n9x.  Now,  Eqn.  (4.3)  is  consistent  with  the  spatial  frequency  com¬ 
ponent  in  [5].  Then,  the  nth  phase  delay  is 

u)0r'n  =  27 rd[n] 

=  -2n(ntfx).  (4.4) 

Proceeding  through  the  radar  model  framework  shows  that  the  return  from  the  rnth 
pulse  at  the  nth  element  is 


%nm 


ate-i2'K’&\n\ei2mn^ 

atej2n(n$x)  ej2itmu> 


(4.5) 


The  spatial  steering  vector  c (9,  0)  must  also  be  verified  using  uniform  element 
spacing.  Inserting  uniform  element  spacing  into  c (0,0)  gives 

c (0,0)  =  e~j2nm  ... 

—  jgjMOD*  ej2-K{l)dx  _  _  .  ej27r(Af-l)^jT  ^  g) 

Equation  (4.6)  is  equivalent  to  the  azimuth  spatial  steering  vector  a(i9x)  defined  in  [5]. 
Changing  c(d,  0)  to  a ($x)  brings  about  many  changes  that  ripple  through  the  radar 
model.  The  3D  steering  vector  defined  in  Chapter  III  as  v(o>,  9 , 0)  =  b(ca)  <g)  c(d,  0) 
can  now  be  written  as 

v(a>,  0)  =  b(u>)  ®  a(i?x).  (4.7) 

The  3D  steering  vector  defined  in  [5]  is 

v(d>,  9 ,  0)  =  e(9z)  ®  b(u>)  ®  a($x)  (4.8) 


where  the  elevation  steering  vector  is 

e(9z)  =  [1  ej2n*z  . . .  e,'2w(p"1)*‘]  . 


(4.9) 
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The  elevation  steering  vector  becomes  e($2)=l  when  P  —  1.  Therefore,  Eqn.  (4.8) 
reduces  to  Eqn.  (4.7)  when  the  uniformly  spaced  array  is  linear.  The  space-time 
snapshot  Xt  has  the  exact  same  change  as  the  3D  steering  vector  since  the  3D  steering 
vector  is  contained  within  the  space-time  snapshot. 

The  next  major  change  to  be  verified  is  the  array  factor  previously  defined  as 
W(0,(f))  =  c(9,(ft)Hc (0',<f/).  The  array  factor  now  becomes  W{4>)  =  a('^)//a('d(.). 
The  array  factor  defined  in  [5] 

W(0,(/>)  =  [e(i?2)  ®  a(-dx)]Hl  (4.10) 

is  for  a  mainbeam  at  boresight.  Equation  (4.10)  is  equivalent  to 

W(0,(/) )  =  [e(d2)  (8)  a(tfx)]H[e(i9'z)  ®  a(d(.)]  (4.11) 

when  the  mainbeam  is  at  boresight.  Equation  (4.11)  simply  allows  the  mainbeam  to 
be  steered  anywhere  and  becomes 

W(4>)  =  a(0x)*a  m  (4.12) 

when  the  array  is  linear.  The  equality  of  Eqns.  (4.10)  and  (4.11)  is  demonstrated  by 
inserting  0  =  0  and  0  =  0  into  '0'x  and  i)'zl  yielding 

Q/  _  dx  cos(0)  sin(0)  _ 

^  X  \  ^ 

Ao 

and 

dz  sin(O) 

U  Z  \ 

^ o 
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Then,  the  elevation  steering  vector  becomes 


e«)  = 


1  e 


J  271-1?'.. 


=  [1 
=  [1  1 


Di27T(0) 


1] 


0j2ir(P-l)d'z 


aj2n(P- 


i)(or 


and  the  azimuth  steering  vector  becomes 


a(^)  = 


1  e 


j2ir&'x 


=  [i 


JMO) 


oj2ir(N—l)'0lx 


J2tt(N- 


1)(0)‘ 


=  [i  i  •••  i] 


(4.13) 


(4.14) 


Since  the  element  pattern  does  not  change  between  each  radar  model  and  the  array 
factors  were  just  shown  to  be  identical  when  the  array  is  linear,  the  antenna  array 
pattern  is  also  equivalent  between  the  two  models. 

The  jammer  and  clutter  models  are  the  final  portion  of  the  radar  model  to 
be  verified.  The  noise  model  is  identical  for  both  uniform  and  non-uniform  element 
spacing  as  described  in  Chapter  III.  The  barrage  noise  jammer  space-time  snapshot 
Xj  =  OLj  (g)  c(9,  (ft)  is  now 

Xj  =  cxj®  a(0x).  (4.15) 

The  barrage  noise  jammer  space-time  snapshot  for  uniform  element  spacing  Xj  = 
e(d2)  <g>  OLj  ®  a(i?x)  reduces  to  Eqn.  (4.15)  when  the  array  is  linear.  The  jammer 
covariance  matrix  changes  to 

R j  =  cr2^jIM  ®  a(^)a^(dx)  (4.16) 

based  on  the  change  in  \y  This  jammer  covariance  matrix  is  equivalent  to  the  one 
defined  in  [13]. 
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The  change  in  the  clutter  model  begins  with  the  ikth  clutter  patch  spatial  fre¬ 


quency, 


k(0i, 

<f>k)  • 

D  [n, : 

A  o 

d  x[n 

cos  9 

i  sin  (j)k 

A0 

—rid; 

c  cos  6 

'i  sin  (pk 

(4.17) 


The  spatial  frequency  definition  of  Eqn.  (4.3)  can  be  substituted  into  Eqn.  (4.17) 
to  give  /dlk [n]  =  —n9x.  The  azimuth  steering  vector  a (i9x)  can  be  substituted  for 
c(@i,  4>k)  in  the  clutter  space-time  snapshot  defined  in  Chapter  III  to  give 

Nr~  1  Nc-1 

Xc  =  £  aikh(u;ik)  <g>  a(0x).  (4.18) 

i= 0  k= 0 

Equation  (4.18)  is  consistent  with  the  clutter  space-time  snapshot  defined  in  [5, 13] 
when  the  array  is  linear.  The  clutter  covariance  matrix 

Nr~  1  Nc-1 

Rc  =  £{XcX f}  =  o-2  ^  £ikHuik)bH (Qik)  ®  c(6'i,0fc)c(6»i,0fc)// 

i= 0  k= 0 
Nr— 1  Nc- 1 

=  or2  tikh(uik)bH (tiik)  ®  a(tfx)aH(dx)  (4.19) 

i=0  fc=0 

is  also  consistent  with  the  clutter  covariance  matrix  defined  in  [5, 13]  as 
Nr—  1  Nc-1 

Rc  =  0-2  X]  ^e(d.)eH(d2)  (8)  b(ujik)bH (uik)  <g>  (tfx)  (4.20) 

i= 0  k= 0 

when  the  array  is  linear. 


4-2  Planar  Array 

The  steering  vectors,  space-time  snapshots,  and  covariance  matrices  are  struc¬ 
tured  differently  in  the  radar  model  presented  in  Chapter  III  than  the  model  in  [5] 
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Table  4.1:  This  table  contains  the  parameters  used  for  the  planar  array 
simulations.  The  simulations  are  performed  to  verify  the  new  model  in 
Chapter  III  is  consistent  with  the  original  model  in  [5]. 


Variable 

Value 

Variable 

Value 

M 

32  pulses 

R 

66  km 

N 

36 

Nc 

NM 

Azimuth  Channels 

6 

P 

1 

Elevation  Channels 

6 

7 

-3  dB 

fo 

1240  MHz 

Array  Transmit  Gain 

22  dB 

dx 

A/2 

Element  Pattern 

Cosine 

dz 

A/2 

Element  Gain 

4  dB 

fr 

1984  Hz 

Element  Backlobe  Level 

-30  dB 

Time  Delay 

0.8  /is 

Transmit  Taper 

Uniform 

Pt 

200  kW 

System  Losses  Ls 

3  dB 

B 

800  kHz 

Target  </> 

0° 

Fn  (Noise  Figure) 

3  dB 

Target  9 

0° 

ha  (aircraft  altitude) 
va  (aircraft  velocity) 

3073  m 

120  m/s 

Target  u 

0.25 

when  the  array  is  planar.  Therefore,  computer  simulation  is  done  to  verify  that  both 
models  produce  identical  results  when  the  array  is  planar  and  has  uniform  element 
spacing.  The  array  used  for  model  verification  is  a  6  x  6  uniformly  spaced  array 
containing  36  elements.  Performance  metrics  include  output  SINR  versus  normalized 
Doppler  frequency  plots  and  detection  probability  versus  input  SINR  per  element  per 
pulse  plots.  All  simulation  parameters  in  this  section  are  given  in  Table  4.1. 


4-2.1  Antenna  Array  Pattern.  Prior  to  simulating  the  performance  metrics, 
the  array  pattern  is  plotted  to  ensure  it  is  consistent  between  both  models.  The  array 
pattern  is  converted  from  radar  coordinates  to  cartesian  coordinates  using  each  model 
and  is  plotted  in  Fig.  4.1.  The  array  center  is  located  at  the  origin  (x  —  0,  y  —  0, 
z  =  0).  As  expected,  most  energy  is  focused  straight  out  of  the  array  at  boresight 
with  symmetric  sidelobes  and  nulls.  The  interference  models  are  examined  next  since 
they  are  directly  affected  by  the  array  pattern. 
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Figure  4.1:  Normalized  array  pattern  in  cartesian  coordinates  for  a  uniformly  spaced 
6x6  element  array  using  the  (a)  new  Chapter  Ill  model  and  (b)  the  original  model 
in  [5].  The  array  center  is  located  at  the  cartesian  coordinate  system  origin  in  this 
plot. 

4-2.2  Covariance  Matrices.  Once  the  array  pattern  is  verified,  the  co- 
variance  matrices  are  checked  for  consistency.  The  rank  of  each  covariance  matrix 
is  examined  by  plotting  the  sorted  eigenvalues  of  each  matrix.  The  rank  of  a  ma¬ 
trix  equals  the  number  of  nonzero  eigenvalues.  Since  the  noise  covariance  matrix 
is  identical  in  both  models  as  discussed  in  Chapter  III,  it  is  not  examined.  The 
noise  covariance  matrix  is  simply  a  scaled  identity  matrix  whose  size  is  determined 
by  the  number  of  array  elements  and  pulses.  The  jammer  and  clutter  matrices,  how¬ 
ever,  are  structured  differently  for  each  model  and  must  be  compared.  The  clutter 
matrix  rank  according  to  Matlab®  is  109  for  both  radar  models  and  its  sorted  eigen¬ 
value  plot  is  shown  in  Fig.  4.2.  The  Brennan  clutter  rank  prediction  is  148  for  each 
clutter  matrix.  The  Brennan  rank  prediction  is  defined  for  a  planar  array  in  [3]  as 
rank(Rc)  ps  [N  +  (M  —  l)/ij  min(P,  Nr  +  1)  where  N  is  the  number  of  azimuth  chan¬ 
nels,  P  is  the  number  of  elevation  channels,  Nr  is  the  number  of  ambiguous  range 
rings,  and  (3  was  discussed  in  Chapter  II.  A  dashed  line  is  drawn  at  the  Brennan 
rank  prediction  to  show  where  the  matrix  eigenvalues  should  go  to  zero.  The  jammer 
matrix  rank  is  32  for  both  models  and  its  sorted  eigenvalue  plot  is  shown  in  Fig.  4.3. 
Both  models  produce  very  similar  sorted  eigenvalue  plots  for  each  matrix,  verifying 
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Figure  4.2:  Clutter  matrix  sorted  eigenvalue  magnitudes  for  a  6x6  array 

with  (3  —  1.  The  clutter  matrix  from  each  model  is  plotted  on  the  same  plot 
to  illustrate  they  are  nearly  identical.  The  Brennan  rank  prediction  shown  by 
the  dashed  line  is  148  while  Matlab®  gives  the  matrix  rank  as  109  in  both 
cases. 

that  the  covariance  matrices  in  each  model  are  the  same.  Although  the  eigenvalue 
magnitudes  of  each  model  in  Fig.  4.3  are  not  identical,  the  important  part  of  the  this 
plot  is  that  the  eigenvalue  magnitudes  drop  at  the  same  point. 


4-2.3  STAP  Algorithms.  The  matched  filter,  adaptive  matched  filter,  JDL, 
and  FTS  STAP  methods  must  also  be  verified  since  they  are  used  in  simulations  in  the 
following  chapter.  Output  SINR  and  detection  probability  plots  are  used  to  compare 
each  STAP  method  for  both  models.  All  figures  shown  overlay  the  simulation  results 
from  each  model.  The  matched  filter  response  is  verified  first  as  shown  in  Fig.  4.4.  The 
matched  filter  output  SINR  is  plotted  as  an  ideal  performance  metric.  Each  STAP 
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Figure  4.3:  Jammer  matrix  eigenvalue  magnitudes  for  a  6x6  array.  Each 

matrix  produces  a  very  similar  sorted  eigenvalue  magnitude  plot.  The  matrix 
rank  is  given  as  32  by  Matlab®  in  each  case. 

method  attempts  to  achieve  the  same  performance,  with  reduced  sample  support,  as 
the  matched  filter.  The  two  models  produce  the  same  results  for  the  matched  filter 
output  SINR.  Output  SINR  peaks  at  about  31  dB.  Figure  4.5  shows  the  output  SINR 
versus  normalized  Doppler  frequency  plot  using  the  JDL  method.  The  JDL  algorithm 
results  are  plotted  using  12  x  4  x  6  =  288  and  4  x  3  x  6  =  72  degrees  of  freedom 
in  Doppler,  azimuth,  and  elevation,  respectively.  As  expected,  SINR  degrades  as  the 
number  of  degrees  of  freedom  decreases.  Figure  4.6  is  another  plot  verifying  the  JDL 
algorithms  from  each  radar  model  produce  the  same  results.  Detection  probability 
is  plotted  versus  input  SINR  per  element  per  pulse  for  the  adaptive  matched  Liter 
and  the  JDL  method  using  1000  realizations  of  x  under  estimated  covariance.  The 
false  alarm  probability  Pfa  is  0.01.  The  false  alarm  probability  means  a  target  will  be 
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Figure  4.4:  Output  SINR  plot  using  the  matched  filter  under  known 
clutter  covariance  for  a  6x6  array. 


detected  one  out  of  every  100  times  when  there  really  is  no  target  present.  Again,  the 
adaptive  matched  filter  is  an  ideal  performance  metric  to  evaluate  JDL  performance. 
JDL  performs  slightly  worse  than  the  adaptive  matched  filter,  as  expected.  The  JDL 
algorithm  requires  about  2  dB  more  input  SINR  than  the  adaptive  matched  filter  to 
obtain  the  same  detection  probability.  The  simulation  of  each  model  was  performed 
separately,  producing  unique  realizations  of  X-  Therefore,  covariance  matrix  estimates 
for  each  model  are  not  based  on  identical  realizations  of  Xi  resulting  in  the  small 
variations  seen  in  detection  probability.  However,  the  detection  probability  results 
are  close  enough  for  model  verification. 

The  other  STAP  method  to  be  verified  is  FTS.  Again,  output  SINR  and  detec¬ 
tion  probability  plots  are  used  to  check  that  the  FTS  algorithm  modified  in  Chap- 
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Figure  4.5:  Output  SINR  plot  of  the  JDL  method  under  known  clutter 

covariance  for  a  6x6  array.  The  JDL  algorithm  in  (a)  uses  12  x  4  x  6  =  288 
degrees  of  freedom  and  (b)  uses  4  x  3  x  6  =  72  degrees  of  freedom  in  Doppler, 
azimuth,  and  elevation.  Output  SINR  degrades  as  the  number  of  degrees  of 
freedom  decreases. 


ter  III  is  consistent  with  the  FTS  algorithm  in  [5].  The  simulated  FTS  algorithm 
uses  a  Blackman-Harris  window  across  Doppler.  Figure  4.7  shows  the  output  SINR 
versus  normalized  Doppler  frequency  plot  for  the  FTS  method.  FTS  performance  is 
notably  worse  than  the  matched  Liter,  as  expected,  since  FTS  is  not  adaptive  across 
Doppler.  Figure  4.8  shows  the  detection  probability  versus  input  SINR  plot  verifying 
the  FTS  and  adaptive  matched  Liter.  As  expected,  FTS  does  not  perform  as  well 
as  the  adaptive  matched  Liter.  FTS  requires  about  7  dB  more  input  SINR  than  the 
adaptive  matched  Liter  to  obtain  the  same  detection  probability. 


4-3  Summary 

The  new  radar  model  presented  in  Chapter  III  is  consistent  with  the  model 
in  [6, 13]  for  the  uniformly  spaced  linear  array  and  the  model  in  [5]  for  the  uniformly 
spaced  planar  array.  This  model  can  now  used  with  conhdence  to  simulate  non- 
uniformly  spaced  arrays.  The  next  chapter  uses  the  model  to  simulate  circular  arrays 
and  arrays  with  failed  elements. 
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Detection  Probability 


Input  SINR  Per  Element  Per  Pulse  [dB] 

Figure  4.6:  Detection  probability  for  the  JDL  method  and  adaptive 

matched  filter  under  estimated  clutter  covariance  for  a  6  x  6  array  using 
1000  realizations  of  Pfa  =  0.01,  and  u  =  0.25.  The  JDL  algorithm 
here  uses  4  x  3  x  6  =  72  degrees  of  freedom  in  Doppler,  azimuth,  and 
elevation.  The  JDL  algorithm  performs  very  well,  only  slightly  worse 
than  the  adaptive  matched  filter. 
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Output  SINR  [dB] 
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Figure  4.7:  Output  SINR  plot  of  the  FTS  method  under  known  clutter 
covariance  for  a  6x6  array.  A  Blackman- Harris  window  is  applied  across 
Doppler  in  the  FTS  algorithm  here  to  compensate  for  its  inability  to 
adapt  across  Doppler. 
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Detection  Probability 


Input  SINR  Per  Element  Per  Pulse  [dB] 

Figure  4.8:  Detection  probability  for  adaptive  matched  filter  and 

FTS  method  with  a  Blackman-Harris  window  across  Doppler  under  es¬ 
timated  clutter  covariance  for  a  6x6  array  using  1000  realizations  of 
Pfa  =  0.01,  and  Co  =  0.25.  FTS  needs  about  7  dB  more  input  SINR 
per  element  per  pulse  to  obtain  the  same  detection  probability  as  the 
adaptive  matched  filter. 
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V.  Non-Uniform  Array  Performance  Evaluation 


This  chapter  uses  the  model  presented  in  Chapter  III  to  examine  the  performance 
of  certain  non-uniformly  spaced  arrays  using  STAP  techniques.  There  are  in¬ 
finitely  many  array  configurations,  so  the  circular  array  is  one  array  configuration 
explored  since  it  has  been  researched  in  the  literature.  The  performance  of  each  array 
is  evaluated  based  on  SINR  loss  plots  and  detection  probability  plots  for  the  FTS 
algorithm,  JDL  algorithm,  and  matched  filter.  SINR  loss  is  plotted  instead  of  output 
SINR  since  the  simulated  arrays  have  different  numbers  of  elements.  Output  SINR  in¬ 
creases  when  more  elements  are  in  the  array.  SINR  loss  is  a  measure  of  actual  output 
SINR  to  maximum  output  SNR.  Maximum  output  SNR  refers  to  the  thermal  noise 
only  case.  Therefore,  SINR  loss  characterizes  array  efficiency  regardless  of  the  number 
of  elements  in  the  array.  The  antenna  pattern  of  each  array  is  also  plotted  to  see  the 
radiated  energy  and  examine  any  sidelobes  present.  Power  spectral  density  (PSD) 
plots  are  provided  to  illustrate  the  impact  of  non-uniformly  spaced  arrays  on  signal 
detection  in  the  interference  environment.  Finally,  the  susceptibility  of  each  array 
to  element  failure  is  examined.  Output  SINR  is  plotted  versus  normalized  Doppler 
frequency  to  study  array  failure  effects.  Certain  non-uniformly  spaced  array  config¬ 
urations  may  potentially  incur  less  degradation  than  uniformly  spaced  planar  arrays 
when  elements  fail.  A  10%  and  20%  element  failure  rate  is  applied  to  each  array  and 
compared  to  array  performance  without  element  failure.  Each  non-uniformly  spaced 
array  is  compared  to  a  6  x  6  uniform  planar  array  with  half  wavelength  inter-element 
spacing.  The  6x6  array  and  non-uniformly  spaced  arrays  have  equal  physical  dimen¬ 
sions,  60.484  centimeters  wide  and  60.484  centimeters  high.  Comparisons  are  done 
to  see  if  the  circular  arrays  perform  as  well  as  the  uniform  6x6  planar  array  having 
the  same  width  and  height.  Non-uniform  element  spacing  may  allow  the  use  of  fewer 
elements  in  the  array  without  sacrificing  signal  detection  performance.  All  simula¬ 
tion  parameters  in  this  chapter  are  the  same  as  the  parameters  given  in  Table  4.1  in 
Chapter  IV  unless  otherwise  specified. 
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5.1  Circular  Array  -  24  Elements 

5.1.1  Configuration.  The  first  non- uniformly  spaced  array  being  examined 
is  a  circular  array  containing  24  elements.  The  array  is  60.484  centimeters  wide  and 
60.484  centimeters  high.  The  distance  between  each  element  is  about  7.87  centime¬ 
ters.  Ideally,  the  elements  should  be  spaced  by  A/2  =  12.10  centimeters  to  prevent 
grating  lobes.  The  element  configuration  of  the  24  element  circular  array  is  shown 
in  Fig.  5.1.  The  platform  velocity  vector  is  along  the  positive  x-axis  since  the  ar¬ 
ray  is  side-looking.  The  array  extends  into  the  negative  x-axis  and  negative  z-axis. 
The  reference  element  in  this  circular  array  is  located  at  the  coordinate  (x  =  0  m, 
z  =  —0.3024  m).  The  reference  element  is  located  furthest  to  the  left,  or  closest 
to  the  aircraft  nose.  The  24  elements  are  arranged  such  that  there  are  13  elevation 
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Figure  5.1:  Element  positions  of  a  side-looking  circular  array  contain¬ 
ing  24  elements.  The  platform  velocity  vector  is  along  the  positive  x-axis, 
hence  it  is  a  sidelooking  array. 
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channels  and  13  azimuth  channels.  The  6x6  array  has  6  elevation  channels  and  6 
azimuth  channels.  The  circular  array  has  better  angular  resolution  than  the  6x6 
array  since  the  number  of  azimuth/elevation  channels  determines  how  well  a  target 
can  be  resolved  in  azimuth/elevation.  For  example,  an  array  with  only  a  few  eleva¬ 
tion  channels  might  resolve  a  target  within  ±15°  of  its  actual  elevation  angle  while  an 
antenna  with  more  elevation  channels  might  be  able  to  resolve  a  target  within  ±5°  of 
its  actual  elevation  angle. 

5.1.2  Array  Pattern.  The  array  pattern  for  the  circular  array  is  shown  in 
Fig.  5.2.  The  array  pattern  is  normalized  and  presented  in  cartesian  coordinates  since 


x-axis  [dB]  60  y-axis  [dB] 


Figure  5.2:  Normalized  3D  array  pattern  converted  from  radar  coordi¬ 

nates  to  cartesian  coordinates  for  the  24  element  circular  array.  The  array 
center  is  located  at  the  origin  (x  —  0,  y  —  0,  z  —  0).  The  element  backlobe 
scaling  factor  is  be  =  —  30  dB. 
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cartesian  coordinates  give  a  more  intuitive  view  of  the  pattern  than  radar  coordinates. 
The  sidelobes  and  nulls  are  circular  in  shape  and  wrap  around  the  mainbeam.  The 
sidelobes  are  about  9  dB  below  the  mainbeam.  The  large  sidelobes  are  undesirable 
because  they  illuminate  clutter  and  other  unwanted  interference,  making  target  detec¬ 
tion  more  difficult.  The  small  number  of  nulls  is  another  disadvantage  for  the  circular 
array  because  only  a  few  azimuth  and  elevation  angles  can  be  nulled  out  to  reduce 
interference,  whereas  the  6x6  uniform  array  in  Fig.  4.1  has  many  more  nulls. 

5.1.3  Clutter  Covariance  Matrix.  The  clutter  matrix  sorted  eigenvalue  mag¬ 
nitudes  for  the  circular  array  and  uniform  planar  array  are  plotted  to  compare  the 
impact  of  each  array  configuration  on  clutter  rank.  The  matrix  rank  depends  on 


Figure  5.3:  Clutter  matrix  sorted  eigenvalue  magnitudes  for  the  24 

element  circular  array  and  the  6x6  planar  array  with  (3  —  1.  The  clutter 
rank  is  similar  for  both  arrays,  meaning  that  the  amount  of  interference 
due  to  clutter  is  only  slightly  different  for  each  array. 
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the  number  of  nonzero  eigenvalues.  Figure  5.3  shows  the  clutter  rank  comparison 
between  the  circular  array  and  the  uniform  planar  array.  Clutter  rank  is  not  signifi¬ 
cantly  higher  for  the  circular  array  when  compared  to  the  uniform  planar  array  in  this 
case.  Therefore,  the  amount  of  interference  due  to  clutter  is  almost  the  same  for  both 
arrays.  Clutter  is  just  as  difficult  to  suppress  with  the  circular  array  as  the  uniform 
planar  array  since  their  clutter  ranks  are  so  similar.  The  slight  difference  in  clutter 
rank  between  the  two  arrays  can  be  attributed  to  a  few  different  factors.  First,  more 
sidelobes  from  the  circular  array  pattern  illuminate  more  clutter  than  the  uniform 
planar  array.  The  circular  array  also  has  fewer  nulls  than  the  6x6  array,  so  more 
interference  is  nulled  out  by  the  uniform  planar  array.  Another  factor  contributing 
to  the  circular  array’s  larger  clutter  rank  is  the  element  spacing.  The  6x6  array  has 
uniform  inter-element  spacing  every  half  wavelength.  Therefore,  clutter  samples  are 
duplicated  since  each  element  moves  a  half  wavelength  per  PRI.  The  circular  array 
does  not  have  uniform  spacing,  so  new  clutter  samples  are  introduced  each  PRI. 

5.1.4  Power  Spectra.  Figure  5.4  shows  the  signal  match  power  spectra  of 
the  24  element  circular  array  and  the  uniform  planar  array.  The  signal  match  plot 
shows  the  signal  received  without  any  adaptivity.  In  other  words,  the  signal  match 
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Figure  5.4:  Signal  match  power  spectra  for  the  (a)  24  element  circular  array 

using  32  pulses  and  the  (b)  6  x  6  planar  array  using  32  pulses.  There  is  more 
interference  introduced  by  the  large  sidelobes  from  the  circular  array. 
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Figure  5.5:  Minimum  variance  estimator  power  spectra  for  the  (a)  24  element 
circular  array  using  32  pulses  and  the  (b)  6  x  6  planar  array  using  32  pulses. 


power  spectra  illustrates  the  signal  return  from  the  environment  based  purely  on  the 
array  pattern.  The  signal  match  power  spectra  shows  how  interference  such  as  clutter 
and  noise  is  illuminated  by  the  mainbeam  and  sidelobes.  Comparing  the  two  power 
spectra  shows  that  more  interference  is  introduced  into  the  system  when  the  array 
is  circular.  More  interference  results  from  the  larger  sidelobes  in  the  circular  array 
pattern  and  the  non-uniform  element  spacing. 

The  MVE  power  spectra  is  shown  in  Fig.  5.5  for  each  array.  The  MVE  power 
spectra  simply  illustrates  the  clutter  location.  The  output  power  is  not  shown  in  the 
MVE  power  spectra  because  amplitude  information  is  meaningless  for  this  algorithm. 
Amplitude  information  is  meaningless  because  MVE  sacrifices  amplitude  estimation 
for  high  resolution  interference  location.  As  expected,  most  interference  is  located 
boresight  to  the  array  at  0  Doppler  frequency.  All  other  interference  is  along  the 
clutter  ridge. 

5.1.5  Matched  Filter.  SINR  loss  versus  normalized  Doppler  frequency  plots 
are  generated  for  the  circular  array.  The  6x6  planar  array  SINR  loss  is  also  plotted 
along  with  the  circular  array  SINR  loss  to  compare  the  performance  of  each  array. 
SINR  loss  shows  the  degradation  of  SINR  due  to  correlated  interference.  The  ideal 
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Figure  5.6:  SINR  loss  using  the  matched  filter  for  the  24  element 

circular  array  and  the  6x6  uniform  planar  array.  The  SINR  loss  is  very 
similar  between  the  two  arrays. 

SINR  loss  is  0  dB,  however  correlated  sources  such  as  clutter  and  jammers  degrade 
SINR.  Therefore,  the  thermal  noise  only  case  is  a  theoretical  performance  bound  used 
to  measure  the  amount  of  degradation  caused  by  correlated  interference  sources. 

The  first  SINR  loss  plot  is  the  matched  filter  scenario  shown  in  Fig.  5.6.  The 
matched  filter  is  fully  adaptive  and  represents  the  best  performance  possible  for  each 
array.  The  output  SINR  for  the  uniform  planar  array  is  about  2  dB  higher  than  the 
output  SINR  for  the  circular  array  across  all  Doppler  frequencies  since  it  has  12  more 
elements.  SINR  tapers  off  at  the  same  Doppler  frequency  for  each  array,  meaning  the 
minimum  and  maximum  discernable  velocities  are  the  same  for  each  array.  Figure  5.6 
shows  the  24  element  circular  array  and  the  6x6  uniform  array  have  very  similar 
SINR  losses  due  to  interference  when  using  the  matched  filter  STAP  algorithm.  The 
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24  element  circular  array  is  efficient  when  compared  to  the  uniform  6x6  array  because 
it  has  12  fewer  elements,  yet  practically  suppresses  correlated  interference  as  well  as 
the  uniform  planar  array. 

5.1.6  JDL  Algorithm.  The  SINR  loss  plot  in  Fig.  5.7  compares  the  circular 
array  and  uniform  planar  array  using  the  JDL  algorithm.  The  JDL  algorithm  used 
in  this  simulation  has  4  x  3  x  6  =  72  degrees  of  freedom  in  Doppler,  azimuth,  and 
elevation,  respectively.  The  reduced  degrees  of  freedom  degrade  output  SINR,  causing 
higher  SINR  losses,  when  compared  to  the  matched  filter  output  SINR.  The  6x6 


Figure  5.7:  SINR  loss  using  the  JDL  algorithm  for  the  24  element 

circular  array  and  the  6x6  planar  array.  The  JDL  algorithm  implemented 
here  uses  4  x  3  x  6  =  72  degrees  of  freedom  in  Doppler,  azimuth,  and 
elevation,  respectively.  The  6x6  uniform  array  has  about  1  dB  more 
SINR  loss  than  the  circular  array  along  the  clutter  null  and  between  0.3 
and  0.7  normalized  Doppler  frequency. 
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uniform  array  has  slightly  more  SINR  loss  than  the  24  element  circular  array  along  the 
clutter  null  and  between  the  normalized  Doppler  frequencies  0.3  and  0.7.  However, 
the  difference  in  SINR  loss  is  only  about  1  dB  in  those  regions  mentioned.  The  circular 
array  suppresses  correlated  interference  slightly  better  than  the  uniform  planar  array 
when  using  the  JDL  STAP  method.  Therefore,  the  circular  array  suppresses  correlated 
interference  more  efficiently  than  the  uniform  planar  array  using  the  partially  adaptive 
JDL  STAP  algorithm. 

Detection  probability  is  plotted  versus  input  SINR  per  element  per  pulse  in 
Fig.  5.8  for  both  arrays  using  the  JDL  algorithm  and  the  adaptive  matched  filter. 


-65  -60  -55  -50  -45  -40  -35 

Input  SINR  Per  Element  Per  Pulse  [dB] 


Figure  5.8:  Detection  probability  plot  using  the  JDL  algorithm  for  the 
24  element  circular  array  and  the  6x6  planar  array  under  estimated 
clutter  covariance  using  1000  realizations  of  %,  Pfa  =  0.01,  and  to  =  0.5. 
The  JDL  algorithm  in  this  case  uses  4  x  3  x  6  =  72  degrees  of  freedom 
in  Doppler,  azimuth,  and  elevation,  respectively. 
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Again,  the  JDL  algorithm  here  uses  4  x  3  x  6  =  72  degrees  of  freedom  in  Doppler, 
azimuth,  and  elevation,  respectively.  A  normalized  Doppler  frequency  of  0.5  is  used  in 
this  plot.  The  circular  array  has  slightly  less  SINR  loss  than  the  uniform  6x6  array 
at  u>  —  0.5  .  The  plot  shows  the  24  element  circular  array  requires  approximately  1  dB 
more  input  SINR  per  element  per  pulse  to  achieve  the  same  detection  probability  as 
the  6x6  array  for  both  STAP  methods.  Fewer  elements  in  the  circular  array  provide 
less  integration  gain  than  the  uniform  6x6  array,  resulting  in  the  small  difference 
in  detection  probability.  The  adaptive  matched  filter  requires  about  2  dB  less  input 
SINR  per  element  per  pulse  than  the  JDL  algorithm  to  obtain  the  same  detection 
probability  for  both  arrays.  The  JDL  algorithm  provides  good  detection  probability 
for  both  arrays  since  it  performs  so  closely  to  the  adaptive  matched  filter. 

5.1.7  FTS  Algorithm.  Figure  5.9  shows  the  SINR  loss  plot  for  each  array 
using  the  FTS  algorithm.  A  Blackman-Harris  window  is  applied  across  Doppler.  The 
window  attempts  to  improve  performance  since  FTS  cannot  adapt  across  Doppler. 
The  narrow  range  of  discernible  velocities  is  a  result  of  the  FTS  algorithm’s  inability 
to  adapt  across  Doppler.  SINR  loss  for  the  24  element  circular  array  is  about  1  dB 
more  than  the  uniform  6x6  array  SINR  loss  when  using  the  FTS  algorithm.  The 
FTS  method  SINR  loss  is  about  3  dB  more  than  the  JDL  method  SINR  loss  across  its 
discernible  velocity  range.  The  JDL  method  has  a  wider  range  of  discernible  velocities 
than  the  FTS  method  for  both  arrays. 

Detection  probability  is  plotted  in  Fig.  5.10  using  the  FTS  algorithm  and  AMF 
for  both  arrays.  The  24  element  circular  array  requires  about  3  dB  more  input  SINR 
per  element  per  pulse  to  obtain  the  same  detection  probability  as  the  uniform  6x6 
array  when  using  the  FTS  algorithm.  The  FTS  algorithm  requires  about  8  dB  more 
input  SINR  per  element  per  pulse  to  obtain  the  same  detection  probability  as  the 
AMF  for  the  circular  array.  The  6x6  uniform  array  requires  about  6  dB  more 
input  SINR  per  element  per  pulse  when  using  the  FTS  algorithm  to  obtain  the  same 
detection  probability  as  the  AMF.  These  results  indicate  the  24  element  circular  array 
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Figure  5.9:  SINR  loss  using  the  FTS  algorithm  for  the  24  element 

circular  array  and  the  6x6  planar  array.  A  Blackman-Harris  window 
is  applied  across  Doppler  since  FTS  cannot  adapt  across  Doppler.  The 
circular  array  has  about  1  more  dB  of  SINR  loss  than  the  6x6  uniform 
array  across  the  discernible  Doppler  frequency  range. 

suffers  slightly  more  than  the  uniform  planar  array  when  using  the  partially  adaptive 
FTS  STAP  method.  Input  SINR  per  element  per  pulse  must  be  increased  by  roughly 
7  dB  using  the  FTS  algorithm  to  obtain  the  same  detection  probability  as  the  JDL 
algorithm  for  the  circular  array.  The  uniform  6x6  array  requires  about  5  dB  more 
input  SINR  per  element  per  pulse  when  using  the  FTS  algorithm  to  obtain  the  same 
detection  probability  when  using  the  JDL  algorithm.  The  reason  FTS  requires  a 
higher  input  SINR  than  the  other  two  algorithms  is  because  FTS  is  not  adaptable 
across  Doppler.  The  FTS  method  cannot  suppress  as  much  interference  as  the  JDL 
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Figure  5.10:  Detection  probability  plot  using  the  FTS  algorithm  and 
AMF  for  the  24  element  circular  array  and  the  6x6  planar  array  under 
estimated  clutter  covariance  using  1000  realizations  of  x,  Pfa  —  0.01,  and 
u  =  0.5.  The  FTS  algorithm  here  employs  a  Blackman-Harris  window 
across  Doppler. 

method.  Therefore,  signal  strength  must  be  increased  if  the  FTS  algorithm  is  to  be 
used. 


5.1.8  Summary.  The  24  element  circular  array  performance  is  similar  to  the 
6x6  uniform  array.  The  circular  array  clutter  rank  is  slightly  higher  than  the  uniform 
planar  array  clutter  rank  due  to  the  larger  sidelobes  and  non-uniform  spacing.  SINR 
loss  is  very  similar  for  both  arrays  when  using  the  matched  filter.  SINR  loss  is  actually 
higher  for  the  uniform  6x6  array  than  the  24  element  circular  array  when  using  the 
JDL  algorithm.  This  result  indicates  the  circular  array  mitigates  correlated  interfer¬ 
ence  more  efficiently  than  the  uniform  6x6  array  under  the  simulation  parameters. 
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The  24  element  circular  array  requires  about  1  dB  more  input  SINR  per  element  per 
pulse  than  the  uniform  planar  array  to  obtain  the  same  detection  probability  as  the 
uniform  planar  array  when  using  the  JDL  method.  The  circular  array  provides  corre¬ 
lated  interference  mitigation  and  detection  probability  comparable  to  the  uniformly 
spaced  6x6  array  using  the  JDL  and  AMF  STAP  techniques.  The  circular  array 
suffers  higher  losses  than  the  uniform  6x6  array  when  using  the  FTS  technique.  The 
performance  metrics  clearly  demonstrate  the  24  element  circular  array  is  an  efficient 
array  when  compared  to  the  uniform  6x6  array. 

5.2  Circular  Array  -15  Elements 
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Figure  5.11:  Element  positions  of  a  side-looking  circular  array  contain¬ 
ing  15  elements.  The  platform  velocity  vector  is  along  the  positive  x-axis. 
The  spacing  between  each  element  is  about  half  a  centimeter  longer  than 
half  a  wavelength. 
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5.2.1  Configuration.  The  next  array  compared  to  the  6x6  uniform  planar 
array  is  a  15  element  circular  array  shown  in  Fig.  5.11.  Both  arrays  have  the  same 
width  and  height  again.  The  15  element  circular  array  is  chosen  because  the  distance 
between  each  element  is  very  close  to  half  a  wavelength.  In  fact,  the  distance  be¬ 
tween  each  element  is  almost  half  of  a  centimeter  longer  than  half  a  wavelength.  The 
inter-element  spacing  is  chosen  to  prevent  grating  lobes.  The  element  positions  are 
arranged  such  that  there  are  9  different  channels  in  azimuth  and  15  different  channels 
in  elevation.  The  15  element  circular  array  pattern  is  not  shown  because  it  is  prac¬ 
tically  identical  to  the  24  element  circular  array  pattern.  The  sidelobes  are  circular 
and  are  about  9  dB  below  the  mainbeam.  The  clutter  rank  plot  is  also  approximately 
the  same  for  the  two  different  circular  arrays.  The  15  element  circular  array  clutter 
rank  is  only  slightly  higher  than  the  uniform  planar  array  clutter  rank.  Larger  side- 
lobes  and  non-uniformly  spaced  elements  increase  clutter  rank  by  introducing  more 
interference  from  clutter. 

5.2.2  Power  Spectra.  The  signal  match  and  MVE  power  spectra  for  the  15 
element  circular  array  and  the  uniform  planar  array  are  compared  in  Figs.  5.12  and 
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Figure  5.12:  Signal  match  power  spectra  for  the  (a)  15  element  circular  array  using 
32  pulses  and  the  (b)  6  x  6  planar  array  using  32  pulses.  The  larger  sidelobes  on 
the  circular  array  illuminate  more  interference  in  the  environment  than  the  uniform 
planar  array. 
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Figure  5.13:  Minimum  variance  estimator  power  spectra  for  the  (a)  15  element 
circular  array  using  32  pulses  and  the  (b)  6  x  6  planar  array  using  32  pulses. 


5.13.  The  circular  array  introduces  more  interference  than  the  uniform  planar  array 
because  of  its  larger  sidelobes  and  non-uniform  element  spacing.  The  signal  match 
output  power  spectra  for  the  15  element  circular  array  is  below  that  of  the  24  element 
circular  array.  This  result  makes  sense  since  output  SINR  increases  as  the  number  of 
elements  in  the  array  increases.  The  15  element  circular  array  MVE  power  spectra  is 
nearly  identical  to  the  24  element  circular  array  MVE  power  spectra.  The  majority  of 
interference  is  located  at  0  normalized  Doppler  frequency  and  0°  azimuth.  All  other 
interference  is  along  the  clutter  ridge. 

5.2.3  Matched  Filter.  The  matched  filter  SINR  loss  is  plotted  in  Fig.  5.14 
using  the  15  element  circular  array  and  the  uniform  planar  array.  The  circular  ar¬ 
ray  output  SINR  is  about  4  dB  below  the  uniform  planar  array  output  SINR  due  to 
the  different  number  of  elements.  However,  SINR  loss  from  correlated  interference 
is  approximately  the  same  for  the  two  arrays  across  all  Doppler  frequencies.  The  15 
element  circular  array  suppresses  correlated  interference  very  efficiently  when  com¬ 
pared  to  the  uniform  planar  array  using  the  matched  filter  since  both  arrays  have 
approximately  the  same  SINR  losses  yet  the  circular  array  has  21  fewer  elements. 
SINR  loss  ranges  from  0  to  4  dB  across  the  discernible  velocity  range  due  to  Doppler 
filter  straddling  losses. 
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Figure  5.14:  SINR  loss  using  the  matched  filter  for  the  15  element 

circular  array  and  the  6x6  uniform  planar  array.  SINR  losses  are  very 
similar  for  both  arrays. 

5.2.4  JDL  Algorithm.  The  JDL  algorithm  is  used  for  the  SINR  loss  plot 
in  Fig.  5.15.  There  are  4  x  3  x  6  =  72  degrees  of  freedom  in  Doppler,  azimuth,  and 
elevation,  respectively.  The  circular  array  suppresses  correlated  interference  quite 
well  when  compared  to  the  uniform  planar  array  using  the  JDL  algorithm.  The 
uniform  planar  array  loses  about  1  dB  more  SINR  than  the  15  element  circular  array 
around  the  clutter  null  and  between  normalized  Doppler  frequencies  0.3  and  0.7.  The 
15  element  circular  array  suppresses  correlated  interference  better  than  the  uniform 
planar  using  the  partially  adaptive  JDL  technique.  Also,  the  JDL  algorithm  results 
in  roughly  1  dB  more  SINR  loss  than  the  matched  filter  for  both  arrays. 

Detection  probability  is  plotted  in  Fig.  5.16  for  the  15  element  circular  array 
and  uniform  planar  array  using  the  JDL  and  AMF  algorithms.  There  are  4x3x6  = 
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Figure  5.15:  SINR  loss  using  the  JDL  algorithm  for  the  15  element 

circular  array  and  the  6x6  planar  array.  The  JDL  algorithm  implemented 
here  uses  4  x  3  x  6  =  72  degrees  of  freedom  in  Doppler,  azimuth,  and 
elevation,  respectively.  The  uniform  planar  array  has  about  1  dB  more 
SINR  loss  than  the  circular  array  between  0.3  and  0.7  normalized  Doppler 
frequency. 

72  degrees  of  freedom  in  Doppler,  azimuth,  and  elevation,  respectively.  Detection 
probability  is  evaluated  at  a  normalized  Doppler  frequency  of  0.5.  At  0.5  normalized 
Doppler  frequency,  the  uniform  6x6  array  has  slightly  more  SINR  loss  than  the 
circular  array.  The  circular  array  requires  about  4  dB  more  input  SINR  per  element 
per  pulse  to  obtain  the  same  detection  probability  as  the  uniform  6x6  array  when 
using  the  AMF.  The  circular  array  requires  about  3  dB  more  input  SINR  per  element 
per  pulse  to  obtain  the  same  detection  probability  as  the  uniform  planar  array  using 
the  JDL  algorithm.  These  results  indicate  the  circular  array  performs  well  when  using 
the  JDL  STAP  algorithm. 
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Figure  5.16:  Detection  probability  plot  using  the  JDL  algorithm  for 

the  15  element  circular  array  and  the  6x6  planar  array  under  estimated 
clutter  covariance  using  1000  realizations  of  %,  Pfa  =  0.01,  and  Q  =  0.5. 

The  JDL  algorithm  in  this  case  uses  4  x  3  x  6  =  72  degrees  of  freedom 
in  Doppler,  azimuth,  and  elevation,  respectively. 

5.2.5  FTS  Algorithm.  The  FTS  algorithm  is  used  for  the  SINR  loss  plot  in 
Fig.  5.17  for  the  15  element  circular  array  and  the  uniform  planar  array.  A  Blackman- 
Harris  window  is  applied  across  Doppler  again.  SINR  loss  for  the  15  element  circular 
array  is  about  the  same  as  SINR  loss  for  the  24  element  circular  array  when  using  the 
FTS  method.  SINR  loss  for  the  15  element  circular  array  is  about  1  dB  lower  than 
the  uniform  6x6  array.  There  is  about  4  to  6  dB  SINR  lost  across  the  discernible 
velocity  range  for  the  circular  array  due  to  correlated  interference.  As  expected,  the 
discernible  velocity  range  is  quite  narrow  when  using  the  FTS  algorithm  because  FTS 
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Figure  5.17:  SINR  loss  using  the  FTS  algorithm  for  the  15  element 

circular  array  and  the  6x6  planar  array.  A  Blackman-Harris  window  is 
applied  across  Doppler  since  FTS  cannot  adapt  across  Doppler. 

is  not  adaptable  in  Doppler.  The  discernible  velocity  range  is  roughly  0.2  to  0.8 
normalized  Doppler  frequency. 

Detection  probability  is  plotted  in  Fig.  5.18  using  the  FTS  and  adaptive  matched 
filer  algorithms  for  both  arrays.  The  FTS  algorithm  performs  notably  worse  than  the 
adaptive  matched  filter  for  both  arrays.  The  FTS  method  requires  about  7  dB  more 
input  SINR  per  element  per  pulse  to  obtain  the  same  detection  probability  as  the 
AMF  for  the  uniform  6x6  array.  For  the  15  element  circular  array,  the  FTS  method 
requires  about  9  dB  more  input  SINR  per  element  per  pulse  to  obtain  the  same 
detection  probability  as  the  AMF.  The  circular  array  suffers  more  than  the  uniform 
6x6  array  when  using  the  FTS  algorithm. 
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Figure  5.18:  Detection  probability  plot  using  the  FTS  algo¬ 
rithm  for  the  15  element  circular  array  and  the  6x6  planar  array 
under  estimated  clutter  covariance  using  1000  realizations  of  x, 

Pfa  =  0.01,  and  Co  =  0.5.  The  FTS  algorithm  here  employs  a 
Blackman- Harris  window  across  Doppler. 

5.2.6  Summary.  The  15  element  circular  array  performance  is  very  similar 
to  the  24  element  circular  array  performance.  The  clutter  rank  for  the  15  element 
circular  array  is  slightly  higher  than  the  uniform  6x6  array  clutter  rank.  Both  arrays 
have  the  practically  the  same  SINR  loss  when  using  the  matched  filter.  When  using 
the  JDL  algorithm,  SINR  loss  for  the  uniform  6x6  array  is  roughly  1  dB  higher 
than  the  15  element  circular  array  SINR  loss  between  0.3  and  0.7  normalized  Doppler 
frequency.  The  15  element  circular  array  requires  about  3  dB  more  input  SINR  per 
element  per  pulse  to  obtain  the  same  detection  probability  as  the  uniform  6x6  array 
using  the  JDL  algorithm.  Under  the  FTS  method,  the  15  element  circular  array 
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does  not  perform  as  well  as  the  uniform  planar  array.  The  15  element  circular  array 
mitigates  correlated  interference  as  well  as  the  uniform  6x6  array  using  the  matched 
filter  and  JDL  STAP  techniques.  The  15  element  circular  array  also  provides  good 
detection  probability  considering  the  small  number  of  elements.  Therefore,  results 
indicate  the  15  element  circular  array  performs  efficiently  when  compared  to  the 
uniform  6x6  array. 

5.3  Array  Failure 

Array  failure  is  briefly  explored  here  since  array  performance  may  vary  according 
to  array  configuration.  Output  SINR  averaged  over  500  realizations  is  plotted  versus 
normalized  Doppler  frequency  using  the  matched  filter  to  characterize  the  performance 
of  each  array.  The  failed  elements  are  randomly  chosen  for  each  realization.  The 
purpose  here  is  to  determine  the  susceptibility  of  each  array  to  element  failure.  In 
all  simulations,  the  array  gain  is  kept  the  same,  meaning  SNR  is  constant  despite 
failing  elements.  In  reality,  SNR  would  decrease  if  some  elements  failed  to  transmit 
because  signal  power  would  decrease.  However,  a  constant  SNR  isolates  degradations 
caused  by  changes  in  the  array  pattern  such  as  grating  lobes,  a  wider  mainbeam, 
larger  sidelobes,  and  different  sidelobe  and  null  locations.  Different  clutter  patches 
are  illuminated  as  a  result  of  element  failure.  Element  failure  rates  of  approximately 
10%  and  20%  are  used  in  all  simulations.  Actual  failure  rates  depend  on  the  number 
of  elements  in  the  array.  A  more  robust  array  will  not  be  significantly  impacted  when 
some  of  its  elements  fail  to  operate. 

5.3.1  6x6  Array.  The  first  array  examined  is  the  6x6  uniformly  spaced 

array.  Average  output  SINR  is  plotted  in  Fig.  5.19  when  0,  4,  and  7  elements  fail. 
Actual  failure  rates  are  11%  and  19%  when  4  and  7  elements  fail,  respectively.  Output 
SINR  peaks  at  about  30.5  dB  when  no  elements  fail.  Output  SINR  degrades  by  about 
1  dB  across  the  discernible  velocity  range  for  approximately  10%  of  element  failure. 
Output  SINR  degrades  at  a  quicker  rate  around  the  clutter  null  with  element  failure. 
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Figure  5.19:  Average  output  SINR  from  the  matched  filter  is  plotted 
for  the  uniformly  spaced  6x6  array  with  0,  4,  and  7  failed  elements.  Out¬ 
put  SINR  is  averaged  over  500  hundred  realizations  with  the  designated 
number  of  random  elements  failing  during  each  realization  and  SNR  = 

1. 

A  Doppler  frequency  of  0.04  sees  about  a  4  dB  degrade  in  output  SINR  from  the  0% 
to  11%  element  failure  rate  and  about  2.5  dB  from  the  11%  to  19%  failure  rate. 

5.3.2  24  Element  Circular  Array.  The  next  array  examined  is  the  24  ele¬ 
ment  circular  array.  Figure  5.20  shows  average  output  SINR  plotted  when  0,  3,  and 
5  elements  fail.  Actual  failure  rates  are  13%  and  21%  when  3  and  5  elements  fail, 
respectively.  Output  SINR  peaks  at  about  28.5  dB  with  no  elements  failing.  Output 
SINR  degrades  by  about  1.3  dB  across  the  discernible  velocity  range  from  the  0%  to 
13%  element  failure  rate.  From  a  13%  to  21%  element  failure  rate,  the  24  element 
circular  array  output  SINR  degrades  about  1  dB.  At  the  clutter  null,  output  SINR 
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Figure  5.20:  Average  output  SINR  from  the  matched  filter  is  plotted 
for  the  24  element  circular  array  with  0,  3,  and  5  failed  elements,  as 
labelled.  Output  SINR  is  averaged  over  500  hundred  realizations  with 
the  designated  number  of  random  elements  failing  during  each  realization 
and  SNR  =  1. 

degrades  at  the  same  rate  as  the  6x6  array.  The  non-uniform  circular  configura¬ 
tion  introduces  fluctuations  in  peak  output  SINR  across  all  Doppler  frequencies  with 
element  failure. 

5.3.3  15  Element  Circular  Array.  The  last  array  examined  is  the  15  element 

circular  array.  Figure  5.21  shows  average  output  SINR  plotted  when  0,  2,  and  3 
elements  fail.  Actual  failure  rates  are  13%  and  20%  when  3  and  5  elements  fail, 
respectively.  Output  SINR  degradation  is  quite  similar  to  the  24  element  circular 
array.  With  a  0%  element  failure  rate,  output  SINR  peaks  at  about  26.5  dB  for 
the  15  element  circular  array.  Output  SINR  degrades  by  about  1.4  dB  across  the 
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Figure  5.21:  Average  output  SINR  from  the  matched  filter  is  plotted 
for  the  15  element  circular  array  with  0,  2,  and  3  failed  elements,  as 
labelled.  Output  SINR  is  averaged  over  500  hundred  realizations  with 
the  designated  number  of  random  elements  failing  during  each  realization 
and  SNR  =  1. 

discernible  velocity  range  from  the  0%  to  13%  element  failure  rate.  The  15  element 
circular  array  output  SINR  degrades  by  about  0.8  dB  from  the  13%  to  20%  element 
failure  rate.  At  the  clutter  null,  output  SINR  degrades  at  the  same  rate  as  the  6x6 
array  and  24  element  circular  array. 


5.3.4  Summary.  The  array  failure  plots  demonstrate  that  output  SINR 
performance  degrades  by  approximately  the  same  rate  for  the  uniformly  spaced  6x6 
array  and  both  circular  array  configurations.  Output  SINR  degrades  for  all  three 
arrays  by  roughly  1  dB  from  0%  to  10%  and  10%  to  20%  element  failure  rates.  All 
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three  arrays  also  see  higher  losses  in  output  SINR  at  the  clutter  null  clue  to  element 
failure. 

5.4  Summary 

The  15  element  circular  array  and  24  element  circular  array  have  practically 
the  same  SINR  losses  for  all  three  STAP  algorithms  used.  Detection  probabilities  for 
each  STAP  method  differ  by  roughly  2  dB  input  SINR  per  element  per  pulse  for  the 
two  circular  arrays.  The  SINR  loss  performance  plots  demonstrate  that  an  array  with 
fewer  elements  can  suppress  correlated  interference  as  well  as  an  array  with  more 
elements,  depending  on  the  configuration.  JDL  performs  significantly  better  than 
FTS  for  the  circular  arrays  and  the  uniform  planar  array.  The  JDL  algorithm  has  a 
wider  discernible  velocity  range,  fewer  SINR  losses,  and  greater  detection  probability 
than  the  FTS  algorithm  for  all  arrays.  Therefore,  JDL  is  a  valuable  STAP  method 
to  use  with  circular  arrays  since  it  performs  so  well.  The  circular  arrays  mitigate 
correlated  interference  very  efficiently  using  the  matched  filter  and  JDL  algorithms 
since  they  have  fewer  elements  than  the  uniform  planar  array,  yet  suffer  roughly  the 
same  SINR  losses.  The  circular  arrays  also  suffer  the  same  as  the  uniform  6x6 
array  in  output  SINR  performance  with  element  failure.  Therefore,  the  circular  array 
configuration  is  a  considerable  alternative  to  the  uniformly  spaced  array  when  using 
STAP  techniques  for  airborne  radar  systems. 
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VI.  Conclusions 


An  airborne  radar  model  for  non-uniformly  spaced  arrays  has  been  developed 
and  presented  in  Chapter  III.  The  model  is  based  on  previous  work  in  [5,6,13] 
and  is  valid  for  linear  and  planar  arrays  having  arbitrary  element  locations.  As  it 
stands,  a  3D  array  can  be  simulated  using  the  model,  however  shadowing  effects  have 
not  been  included.  Analysis  of  airborne  radar  systems  with  non-uniformly  spaced 
arrays  is  valuable  since  non-uniformly  spaced  arrays  may  potentially  provide  better 
interference  mitigation  and  target  detection  than  uniformly  spaced  arrays.  The  model 
developed  in  this  thesis  provides  the  mathematical  framework  to  evaluate  any  linear 
or  planar  array  configuration  using  the  matched  filter,  AMF,  FTS,  and  JDL  STAP 
techniques. 

6.1  Airborne  Radar  Model  For  Non-Uniformly  Spaced  Arrays 

The  model  development  begins  by  introducing  a  new  element  position  matrix 
D  to  accommodate  arbitrary  element  locations  in  the  array.  This  element  position 
matrix  creates  a  chain  reaction  affecting  many  components  in  the  airborne  radar 
model.  First,  the  time/phase  delay  to  each  element  changes  since  the  elements  can 
move  to  any  location.  Spatial  frequencies  change  based  on  the  time/phase  delay, 
resulting  in  a  new  spatial  steering  vector  c.  The  array  factor  is  restructured  based 
on  the  new  spatial  steering  vector.  The  space-time  snapshot  is  reorganized  since  it 
contains  the  new  spatial  steering  vector.  Finally,  the  jammer  and  clutter  covariance 
matrices  are  restructured  due  to  the  space-time  snapshots.  Two  STAP  techniques, 
JDL  and  FTS,  are  also  updated  so  they  can  be  applied  to  non-uniformly  spaced 
arrays. 

6.2  Model  Verification 

Verification  of  the  airborne  radar  model  for  non-uniformly  spaced  arrays  was 
accomplished  in  Chapter  IV  by  inserting  uniform  element  spacing  into  the  model. 
The  new  model  presented  in  Chapter  III  mathematically  collapses  back  to  the  orig- 
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inal  model  of  [6, 13]  when  the  array  is  linear.  Mathematical  verification  could  not 
be  accomplished  for  the  uniform  planar  array  due  to  the  restructured  steering  vec¬ 
tors,  space-time  snapshots,  and  covariance  matrices.  Therefore,  simulations  were 
performed  to  compare  the  results  produced  from  the  new  model  in  Chapter  III  and 
the  original  model  in  [5].  Simulations  included  array  pattern  plots,  interference  ma¬ 
trix  eigenvalue  magnitude  plots,  output  SINR  versus  normalized  Doppler  frequency 
plots,  and  detection  probability  versus  input  SINR  per  element  per  pulse  plots.  Each 
model  produced  the  same  results,  verifying  their  equivalence  for  the  uniformly  spaced 
planar  array. 

6.3  Model  Application 

The  model  was  used  to  simulate  a  non-uniformly  spaced  array  after  verification. 
While  the  main  focus  of  this  thesis  was  not  to  determine  the  ideal  array  configuration, 
the  model  presented  in  Chapter  III  provides  the  ability  to  evaluate  the  performance  of 
non-uniformly  spaced  arrays  in  airborne  radar.  Two  different  circular  arrays  were  sim¬ 
ulated  using  this  model,  one  containing  15  elements  and  one  containing  24  elements. 
The  performance  of  each  circular  array  was  compared  to  a  uniformly  spaced  6x6 
array.  All  three  arrays  had  the  same  width  and  height.  It  was  found  that  the  circular 
arrays  introduced  more  interference  from  the  environment  than  the  uniform  planar 
array  clue  to  larger  sidelobes,  fewer  nulls,  and  non-uniform  element  spacing.  However, 
the  three  arrays  had  very  similar  SINR  losses  when  using  the  matched  filter  and  JDL 
algorithms.  The  circular  arrays  actually  had  less  SINR  loss  than  the  uniform  planar 
array  when  using  the  JDL  STAP  algorithm.  The  circular  arrays  had  greater  SINR  loss 
than  the  6x6  array  when  using  the  FTS  algorithm.  Detection  probability  is  difficult 
to  compare  between  the  arrays  since  they  contain  different  numbers  of  elements.  The 
6x6  array  has  more  elements,  so  it  has  a  higher  integration  gain  than  the  circular 
arrays.  However,  detection  probability  results  indicate  JDL  performs  well  with  the 
circular  arrays.  The  FTS  method  did  not  perform  as  well  as  the  JDL  method  for  all 
three  arrays.  The  circular  arrays  degraded  in  output  SINR  at  about  the  same  rate  as 


81 


the  uniform  planar  array  for  the  10%  and  20%  element  failure  rates.  All  simulation 
results  demonstrate  there  are  potential  array  configurations  which  are  more  efficient 
than  the  uniform  planar  or  linear  arrays  typically  used  in  airborne  radar  platforms. 
These  potential  array  configurations  could  lead  to  reduced  array  costs  due  to  fewer 
elements  without  sacrificing  detection  capabilities. 

6.4  Future  Research 

As  it  stands,  the  model  in  Chapter  III  is  valid  for  linear  and  planar  arrays  with 
arbitrary  element  spacing.  The  non-uniformly  spaced  array  model  could  be  modified 
to  include  shadowing  effects  in  the  3D  conformal  array  case.  Incoming  radar  waves 
might  not  strike  every  element  in  the  array  if  it  is  conformal,  so  this  must  be  taken 
into  consideration. 

The  simulations  performed  in  Chapter  V  primarily  demonstrated  the  model 
for  non-uniformly  spaced  arrays  is  valid.  There  are  infinitely  many  arrays  that  can 
be  explored,  so  research  should  be  done  to  investigate  different  array  configurations. 
More  work  is  necessary  to  evaluate  performance  of  each  non-uniform  array.  Detection 
probability  and  SINR  loss  should  be  evaluated  for  a  variety  of  scenarios  to  better 
characterize  overall  array  performance.  Many  different  target  azimuth  and  elevation 
angles,  Doppler  frequencies,  and  ranges  could  be  considered.  The  degrees  of  freedom 
could  also  be  varied  when  using  each  STAP  technique. 

The  non-uniformly  spaced  array  model  can  be  extended  to  different  applications, 
including  space-based  radar  as  one  possible  extension.  Multiple  spacecraft  form  an 
array  in  space,  where  each  spacecraft  is  an  element.  The  space-based  radar  model  will 
have  a  few  changes  that  must  be  made.  The  clutter  model  will  change  since  the  radar 
platform  is  at  a  much  higher  altitude.  Atmospheric  attenuation  must  also  be  included 
for  the  space-based  radar  model.  The  earth’s  rotation  will  also  be  another  factor  to 
be  considered.  There  are  many  different  constellations  presented  in  the  literature  that 
can  be  studied. 


82 


Appendix  A.  Kronecker  Product 


f  j  ^  lie  Kronecker  Product  is  defined  in  [4]  as 

Ai\Q  AV2Q  Ai3B  •  ■  •  AimB 

A21B  A22  B  A23B  •  •  •  42mB 

A  <g>  B  = 

AwiB  Aat2B  Aat3B  •  ■  •  Ajvm B 

Every  element  in  A  is  multiplied  by  the  matrix  B.  If  A  is  N  x  M  and  B  is  P  x  Q, 

then  the  resulting  matrix  is  NP  x  MQ. 


83 


Bibliography 


1.  Balanis  C.  Antenna  Theory  (Second  Edition).  New  York,  NY:  John  Wiley  & 
Sons,  Inc.,  1997. 

2.  Brooks  L.  S.  and  Reed  I.  S.  “Equivalence  of  the  Likelihood  Ratio  Processor,  the 
Maximum  Signal-to-Noise  Ratio  Filter,  and  the  Wiener  Filter,”  IEEE  Transac¬ 
tions  on  Aerospace  and  Electronic  Systems,  690-692  (September  1972).  Corre¬ 
spondence. 

3.  Corbell  P.  M.  and  Hale  T.  B.  “3-Dimensional  STAP  Performance  Analysis  Using 
the  Cross-Spectral  Metric.”  Proceedings  of  the  2004  IEEE  Radar  Conference. 
610-615.  Philadelphia,  PA:  IEEE  Aerospace  and  Electronic  Systems  Society, 
April  2004. 

4.  Graham  A.  Kronecker  Products  and  Matrix  Calculus  with  Applications .  Ellis 
Horwood:  Mathematics  and  its  Applications,  Market  Cross  House,  Cooper  Street, 
Chichester,  West  Sussex,  P019  1EB,  England:  Ellis  Horwood  Limited,  1981. 

5.  Hale  T.  B.  Airborne  Radar  Interference  Suppression  Using  Adaptive  Three- 
Dimensional  Techniques.  Ph.D.  dissertation,  School  of  Engineering  and  Man¬ 
agement,  Air  Force  Institute  of  Technology  (AETC),  2950  P  Street,  Bldg  640, 
Wright-Patterson  AFB,  OH  45433-7765,  June  2002.  AFIT/DS/ENG/02-02. 

6.  Jaffer  A.,  Baker  M.,  Ballance  W.,  and  Staub  J.  Adaptive  Space-Time  Process¬ 
ing  Techniques  for  Airborne  Radars.  Contract  F30602-89-D-0028,  Fullerton,  CA 
92634:  Hughes  Aircraft  Company,  July  1991. 

7.  Klemm  R.  Principles  of  Space- Time  Adaptive  Processing,  12.  IEE  Radar,  Sonar, 
Navigation  and  Avionics  Series.  Michael  Faraday  House,  Six  Hills  Way,  Stevenage, 
Herts.  SGI  2AY,  United  Kingdom:  Institution  of  Electrical  Engineers  (IEE),  2002. 
ISBN  0852961723. 

8.  Nathanson  F.  E.  Radar  Design  Principles:  Signal  Processing  and  the  Environment 
(Second  Edition).  McGraw-Hill,  Inc.,  1991. 

9.  Reed  I.  S.,  Mallett  J.,  and  Brennan  L.  “Rapid  Convergence  Rate  in  Adaptive 
Arrays,”  IEEE  Transactions  on  Aerospace  and  Electronic  Systems,  AES- 10,  No. 
6: 853-863  (November  1974). 

10.  Skolnik  M.  I.  Introduction  to  Radar  Systems  (Third  Edition).  1221  Avenue  of  the 
Americas,  New  York,  NY  10020:  McGraw-Hill,  Inc.,  2001.  ISBN:  0072881380. 

11.  Stutzman  W.  L.  and  Thiele  G.  A.  Antenna  Theory  and  Design  (Second  Edition). 
New  York,  NY:  John  Wiley  &  Sons,  Inc.,  1998. 

12.  Wang  H.  and  Cai  L.  “On  Adaptive  Spatial- Temporal  Processing  for  Airborne 
Surveillance  Radar  Systems,”  IEEE  Transactions  on  Aerospace  and  Electronic 
Systems,  5’0(3):66O-669  (July  1994). 


84 


13.  Ward  J.  Space-Time  Adaptive  Processing  for  Airborne  Radar.  Contract  F19628- 
95-C-0002,  Lexington,  Massachusetts:  Lincoln  Laboratory,  Massachusetts  Insti¬ 
tute  of  Technology,  December  1994. 

14.  Zatman  M.  “Circular  Array  STAP.”  Proceedings  of  the  1999  IEEE  Radar  Con¬ 
ference.  108-112.  April  1999. 

15.  Zatman  M.  “Circular  Array  STAP,”  IEEE  Transactions  on  Aerospace  and  Elec¬ 
tronic  Systems,  36 ( 2):510— 516  (April  2000). 


85 


